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Abstract 


A  computerized  simulation  model  was  developed  as  a  flexible  tool  to 
aid  in  the  design*  deployment,  and  employment  decisions  regarding  High 
Energy  Laser  (HEL)  weapons.  Only  the  "space  defensive"  and  "entisatellite" 
missions  are  considered.  The  simulation  models  the  physical  processes  of 
laser  propagation  and  laser  effects,  and  it  includes  a  "first-order" 
optimization  process  for  selection  of  weapons  and  weapon  paths  of  maximum 
efficiency.  It  also  includes  a  "many-on-many"  engagement  model  to  test  the 
effectiveness  of  the  weapons  selected. 

Target  satellites  include  a  set  of  four  hypothetical  systems  represent¬ 
ing  "real  world"  possibilities.  Weapon  choices  may  be  based  on  the  surface 
of  the  earth,  in  aircraft,  or  in  circular  orbits.  The  user  specifies  weapon 
design  characteristics  as  well  as  an  associated  basing  mode  for  each  type  of 
weapon  to  be  deployed  by  the  optimization  model. 

In  addition  to  selecting  the  targets  and  defining  the  types  of  weapons 
to  be  used,  the  user  defines  the  time  allowed  to  accomplish  a  defined  miss¬ 
ion,  and  the  percentage  of  each  system  of  target  satellites  which  must  be  ne¬ 
gated.  Rudimentary  tactics  for  weapon  employment  are  also  allowed. 

Validation,  verification,  and  experimental  designs  for  model  use  are 
discussed.  R* commendations  for  model  expansion  are  given. 
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I.  INTRODUCTION 


This  thesis  is  an  operations  research  oriented  look  at  part  of 


the  world  of  High  Energy  Laser  (HEL)  weapons  applications.  In  this 


"light,"  the  investigation  of  aspects  such  as  cost  effectiveness,  cost 


benefits,  weapon  trade-offs,  the  strategy  of  deployment,  the  tactics 


of  employment,  mission  accomplishment  and  various  system  optimizations 


are  of  primary  interest. 


This  chapter  includes  sections  which  discuss  the  background 


associated  with  defining  the  need  for  investigation  of  these  things. 


the  general  problem  area  to  which  this  thesis  is  addressed  and  its 


specific  objective,  the  scope  and  limitations  of  the  analysis  and  the 


model,  the  assumptions  made  to  make  the  problem  tractable,  the  approach¬ 


es  to  solving  the  problem  and  building  the  model,  and  the  sequence  of 


chapters  which  follow. 


Background 


Research  and  development  concerning  HEL  devices,  as  would  be 


expected  for  an  infant  technology  with  such  long  range  promise,  has 


experienced  major  funding  and  directional  changes  over  the  past  ten 


years.  The  major  thrust  of  Air  Force  effort  in  the  seventies,  (to 


demonstrate  the  feasibility  of  acquiring,  identifying,  pointing,  track¬ 


ing,  and  disabling  an  aerodynamic  object  from  an  aircraft)  has  apparently 
(Ref  1:32-66)  given  way  to  something  new.  The  "space  defensive"  mission, 


with  a  real  probability  of  a  proof-of- concept  demonstration  being  the 


first  U.S.  HEL  battle  station,  has  apparently  moved  to  a  high  priority 


position  in  the  competition  for  R&D  HEL  resources. 
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The  perception  of  the  level  of  Soviet  effort  in  this  area,  along 
with  their  intentions  and  capabilities,  has  combined  with  a  realization  of 
U.S.  strategic  dependence  on  our  space  systems  to  warrant  this  shift  in 
priorities  within  the  Department  of  Defense  (Ref  2:1-5;  3:317-323). 

As  more  resources  are  expended  in  the  effort  to  field  a  viable  sys¬ 
tem  to  perform  the  space  defensive  mission,  the  decisions  required  in  the 
system  research,  development,  test  and  acquisition  process  become  more 
crucial  (Ref  4:98-100).  To  field  the  right  system  at  the  right  time  in  the 
right  place  and  for  the  right  cost  requires  that  the  best  decisions  possible 
be  made  from  basic  research  through  systems  acquisition.  Toward  this  end, 
a  flexible  tool  in  the  form  of  a  systems  model  has  been  developed. 

Problem 

The  nature  of  the  research  and  development  process  leading  to  proto¬ 
type  hardware  or  proof-of-concept  demonstration  of  HEL  devices  is  certainly 
no  less  risky  than  that  for  other  relatively  new  technologies.  In  the  R&D 
process  for  HEL  systems,  decisions  regarding  the  expenditure  of  resources 
for  basic  research,  laboratory  hardware,  development  of  peripneral  and 
support  equipment,  production  of  weapon  prototypes,  and  the  full  product¬ 
ion  buy  must  be  made  with  some  degree  of  knowledge  concerning  the  probable 
payoffs.  Knowledge  of  probable  payoffs  in  terms  of  beam  power,  beam  size, 
propagation  characteristics,  device,  size  and  weight,  power  requirements, 
etc.  are  certainly  being  weighed  in  the  decision-making  process.  However, 
the  combination  of  these  measurable  weapon  attributes  in  satisfying  a 
specific  military  objective  or  mission  is  perhaps  generally  less  quantified. 

If  the  objective  is  to  develop  and  deploy  an  HEL  weapon  system  for 
satellite  defense,  where  do  we  put  our  money?  We  can  say  that  short 
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wavelengths,  high  efficiencies,  good  reliability,  lightness  in  weight, 
large  optics,  high  pointing  and  tracking  accuracies,  minimum  costs, 
large  energy  storage  capacity  and  the  ability  to  command  the  weapon  are 
all  important  in  the  design  of  HEL  weapons  to  make  up  the  weapon  system. 
But  how  do  we  make  trade-cffs  of  conflicting  attributes:  For  example, 
would  it  be  cheaper  to  expend  resources  to  develop  a  relatively  long 
wavelength  weapon  which  could  be  put  into  several  lower  orbits,  or  to  per¬ 
haps  go  forward  with  research  leading  to  shorter  wavelength  weapons  which 
could  be  put  into  geo- synchronous  orbits?  If  we  are  at  the  prototype 
stage,  and  are  looking  forward  to  a  production  contract,  how  many  weapons 
should  we  buy  and  what  are  the  best  orbits  to  put  them  in?  Would  it  be 
cheaper  to  expend  launch  vehicles  to  place  them  all  at  synchronous  alti¬ 
tude,  or  perhaps  to  place  a  few  more  in  lower  orbits?  Is  there  some 
combination  of  orbits  which  can  satisfy  the  objective  in  some  way  "better" 
than  any  single  orbital  altitude?  And,  going  a  step  beyond  deployment 
options,  what  kinds  of  tactics  are  best  suited  for  satellite  defense? 

The  general  problem  of  interest  here  is  how  best  to  examine  the 
various  attributes  of  HEL  devices  and  deployment  options  as  they  relate 
to  a  specifically  stated  military  objective. 

Thesis  Objective 

The  purpose  of  this  thesis  is  to  properly  define  the  appropri¬ 
ate  attributes  of  merit  and  to  develop  a  simulation  model  which  would 
give  a  decision  maker  a  flexible  tool  in  investigating  the  effects 
of  these  attributes  on  specific  space  defensive  missions. 

Scope  and  Limitations 

Scope.  This  effort  is  limited  to  an  analysis  concerning  satellite 
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targets  and  HEL  devices  based  on  the  surface  of  the  earth,  in  aircraft  or 
in  orbit.  The  purpose  of  the  model  HELBASE  is  to  "most  efficiently"  define 
the  basing  parameters  for  a  constellation  of  weapons  to  meet  the  user-de¬ 
fined  mission  requirements  under  the  constraints  cf  user-supplied  weapons 
characteristics,  allowed  basing  options  and  tactics.  Ground- to-space, 
aircraft-to-space,  and  space-to-space  irradiation  is  allowed.  This  approach 
would  apply  to  investigating  the  characteristics  of  both  satellite  defensive 
systems  or  strictly  anti -satellite  ( AS  AT )  systems. 

The  user  inputs  allowed  are  Target  Type  Selections,  Weapon  Type  Select¬ 
ions,  Mission  Data,  and  Battle  Management  Data. 

The  model's  most  singular  characteristic  is  its  flexibility  as  a  tool 
for  the  operations  analyst.  This  flexibility  as  an  analytical  tool  requires 
a  rather  large  number  of  user  inputs,  and  this  immediately  places  a  number  of 
responsibilities  on  the  user.  Since  virtually  alj_  first-order  variables 
related  to  weapon  system  versus  target  system  engagement  design  are  included 
in  user  inputs,  the  user  must  select  the  control  variables  and  decision 
variables  pertinent  to  the  specific  application  (i.e.  purpose)  the  user 
intends.  In  other  words,  some  thoughtful  experimental  design  is  normally 
necessary  before  employing  HELBASE.  For  instance,  investigation  of  launch 
costs  would  require  controlling  Mission  Data,  Target  Types,  and  probably 
Battle  Management  Data,  and  allowing  the  orbital  selections  of  Weapon  Types 
to  vary  along  a  range  of  interest.  In  this  way,  launch  costs  could  con¬ 
ceivably  be  minimized.  Another  user  may  be  interested  in  evaluating  various 
battle  tactics  or  rules  of  engagement.  This  could  be  done  by  controlling 
Weapon  Types,  Target  Types  and  Mission  Data,  and  then  varying  inputs 
associated  with  Battle  Management  Data. 
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So  this  flexibility  of  HELBASE  will  allow  its  use  in  the  investi¬ 
gation  of  a  wide  range  of  problems  related  to  the  deployment  and  employ¬ 
ment  of  HEL  weapons  against  target  satellites.  This  flexibility  places 
a  responsibility  on  the  user  to  understand  HELBASE  functions  and  methodol¬ 
ogy  before  confidence  can  be  placed  in  the  results.  In  fact,  both  pre- 
and  post-analysis  steps  are  strongly  recommended  for  the  use  of  HELBASE. 
This  will  be  discussed  in  detail  in  Chapter  V. 

Limitations.  It  cannot  be  over-emphasized  that  HELBASE  will  only 
provide  a  first  order  solution  to  the  most  efficient  constellation  of 
allowed  weapon  types  to  perform  the  mission  defined.  The  greatest  degree 
of  approximation  occurs  in  the  modeling  of  laser  propagation  through  the 
atmosphere,  and  in  the  modeling  of  the  effects  these  beams  will  have  on 
the  targets.  This  degree  of  approximation  is  justified  by  the  "first 
order"  accuracy  of  HELBASE,  and  by  the  large  number  of  laser  wavelengths 
allowed. 

Only  circular  orbits  are  allowed  for  the  orbiting  battle  stations 
in  the  general  optimization  algorithm.  Other  than  in  the  placement  of 
HEL  battle  stations  in  the  proximity  of  satellites  to  be  defended,  or 
the  orbit  matching  of  target  satellites  with  highly  elliptical  orbits, 
the  value  of  allowing  elliptical  orbits  for  HEL  battle  stations  is  diffi¬ 
cult  to  substantiate.  In  addition,  the  extra  two  degrees  of  freedom  in 
orbit  selection  (eccentricity  and  argument  of  perigee)  make  the  general 
efficiency  optimization  problem  much  more  formidable  (see  Appendix  B). 
Input  of  specific  elliptical  orbits  for  weapons  would  seem  to  have  value, 
and  is  discussed  in  Appendix  B. 

Obviously,  only  HEL  battle  stations  are  considered,  although 

5 


ir  -  ■mnwnrtmtin-.i  niT—ffrrr  . 


im 


model  expansion  to  include  Charged  Particle  Beam  (CPB)  weapons  would  seem 
worthwhile  (e.g.  trade-off  and  comparisons  between  HEL  and  CPB  weapons  for 
a  specific  mission). 

In  the  same  light,  the  user  will  be  restricted  to  defining  only 
one  weapon  per  battle  station. 

Last,  the  cost  impact  of  launching  the  "most  efficient"  battle 
station  constellation  into  the  required  orbits  is  entirely  ignored,  as 
are  the  fixed  costs  of  using  the  first  weapon  of  a  particular  type  in 
building  the  constellation. 

To  summarize,  HELBASE  is  intended  as  a  first  order  approximation 
to  the  most  efficient  basing  distribution,  and  considers  only  the  satellite 
defense  and  ASAT  missions. 

Assumptions 

Implied  in  the  construction  of  the  HELBASE  optimization  algorithm 
are  the  assumptions  of  zero  launch  and  fixed  production  costs.  Knowledge 
of  this  will  allow  a  judicious  user  to  design  the  model  use  experiment  to 
compensate  if  necessary  for  a  cost  analysis. 

In  consonance  with  a  first-order  approximate  solution,  orbital 
mechanics  assume  point  masses  and  no  perturbations.  Circular  orbits  are 
perfectly  circular,  the  earth  is  a  perfect  sphere  of  radius  6,378.145 
kilometers,  and  so  on. 

Battle  stations  do  not  maneuver  in  orbit  nor  change  location  on 
the  ground  or  in  the  air.  (That  is,  aircraft  carrying  HEL  weapons  are 
assumed  to  "orbit"  a  specific  ground  reference  point). 

There  is  no  "phase  relationship"  between  targets  and  weapons. 
However,  there  ?s  information  available  (this  information  is  fairly  easily 
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obtained)  concerning  phase  relationships  could  be  used  to  improve  the 
deployed  systems  performance.  (See  Chapter  V). 

All  "most  efficient"  orbits  identified  by  HELBASE  are  considered 
feasible  with  available  launch  vehicles  and  sites.  That  is,  user 
specification  of  a  weapon  type  which  involves  placing  10,000  pounds  in 
synchronous  orbit  will  be  considered  feasible  by  HELBASE,  but  would  be 
perhaps  of  small  value  to  the  user. 

The  vulnerabilities  of  targets  are  assumed  to  be  accurate  and 
reliable.  (See  Appendix  A). 

No  thermal  blooming  will  occur  during  atmospheric  propagation 
(See  Chapter  V). 

Last,  and  possibly  most  important,  it  is  assumed  that  the  input 
of  weapon  parameters  by  the  user  implies  that  the  user  can,  at  the  time 
intended,  solve  all  technical  problems  involving  weapon  production, 
required  pointing  and  tracking  accuracies,  etc. 

General  Approach 

In  addition  to  a  description  of  the  general  optimization  function, 
the  approach  to  the  development  of  the  HELBASE  model  will  be  summarized. 

Approach  to  "Most  Efficient"  Weapon  Constellation.  A  functional 
understanding  of  the  process  by  which  HELBASE  defines  the  "most  efficient 
constellation  of  weapons  to  meet  the  mission  requirements  can  be  gained 
by  following  a  simple  case  through  the  optimization  process.  The  range 
of  all  user  inputs  in  the  model  will  be  considered  in  Chapter  II.  Here 
we  will  consider  the  simple  case  of  1  target  type  (say  4  satellites)  and 
1  weapon  type  (say  at  300  nautical  miles)  for  explanatory  purposes. 
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We  first  establish  a  sphere  centered  at  the  earth's  center  whose 
radius  is  the  radius  of  the  earth  plus  300  nautical  miles  (or  555.6 
kilometers).  Our  weapons  must,  by  definition  of  the  weapon  type,  be 
placed  in  orbits  corresponding  to  great  circles  around  this  sphere. 

The  second  step  is  to  calculate  an  average  sighting  density  over  the 
surface  of  this  sphere  by  stochastically  placing  the  target  satellites 
into  their  orbits,  and  determining  whether  they  can  be  "sighted"  from 
each  of  1650  points  distributed  uniformly  over  the  surface  of  the  sphere. 
This  process  is  repeated  many  times  to  establish  the  average  sighting 
density  for  each  point.  Once  these  are  established,  the  great  circle 
which  offers  the  greatest  average  sighting  over  the  circle  is  found. 

A  battle  station  is  then  placed  in  that  great  circle  (circular  orbit). 

A  check  of  whether  this  first  weapon  can  satisfy  the  mission  require¬ 
ments  is  made  by  placing  the  targets  and  this  weapon  stochastically  into 
their  orbits  and  letting  the  weapon  fire  against  the  targets  over  time. 

If,  at  the  end  of  the  Target  System  Negation  Time  (specified  by  the  user), 
all  of  the  mission  requirements  are  met,  the  existing  system  could 
conceivably  meet  the  mission  requirements,  and  no  more  weapons  would  need 
to  be  placed.  If  all  mission  requirements  were  not  met,  we  would  need 
another  weapon.  The  average  sighting  densities  would  then  be  recalculated, 
with  the  effect  of  the  first  weapon  being  accounted  for.  The  entire  pro¬ 
cess  is  repeated  until  the  mission  requirements  can  be  met.  The  result 
is  a  constellation  of  battle  stations  in  300  nautical  mile  orbits  which 
"most  efficiently"  satisfies  the  user's  mission  requirements. 
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Approach  to  Model  Development.  The  modular  concept  of  model  develop¬ 
ment  was  used  to  allow  the  development  of  the  basic  submodels  of  HELBASE 
before  combining  them  into  increasingly  larger  units.  The  basic  development 
steps  are  documented  by  Shannon  (Ref  19:23)  and  others,  and  include: 


System  definition 
Model  formulation 
Data  preparation 
Model  translation 
Verification  and  validation 
Experimental  design 
Experimentation 
Analysis  of  output 
Documentation. 

The  definition  of  the  system  under  consideration  is  discussed  in 
this  chapter,  and  the  logic  flow  reflecting  model  formulation  is  the 
concern  of  Chapter  II.  Data  preparation  involved  some  manipulation  of 
the  data  and  is  explained  in  Appendix  A.  Verification  and  validation 
are  the  subjects  of  Chapter  IV,  and  experimental  design  and  analysis 
of  model  output  are  covered  in  Chapter  V.  Documentation  of  the  HELBASE 
model  is  in  the  form  of  detailed  comments  included  throughout  the  model. 

A  copy  of  the  Fortran  program  is  contained  in  Appendix  C. 

Sequence  of  Presentation 

Chapter  II  (Model  Overview)  discusses  in  detail  the  functions  of 
HELBASE.  Supporting  analyses  for  model  design  and  functions  are  contained 
in  Chapter  III.  Chapter  IV,  Verification  and  Validation,  reviews  the 
procedures  and  checks  which  were  undertaken  to  insure  that  the  HELBASE 
model  will  perform  in  accordance  with  its  stated  purpose  and  within  its 
scope  and  limitations.  Experimental  designs  for  model  use  will  be  dis¬ 
cussed  in  Chapter  V  along  with  some  representative  demonstrations. 
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Conclusions  and  Recommendations,  Chapter  6,  includes  suggestions  for 
model  expansion*  and  these  suggestions  are  discussed  further  in  Appendix 
B. 


II  MODEL  OVERVIEW 


This  chapter  presents  the  general  characteristics  and  capabilities 
of  the  model  HELBASE  by  focusing  on  the  basic  functional  relationships. 
In  addition  to  the  overall  model,  the  user  inputs  and  each  of  the  sub¬ 
models  are  described  in  terms  of  the  general  purpose  and  functions  of 
each.  The  general  underlying  character  of  HELBASE  Is  that  the  output 
weapon  system  will  be  capable  of  meeting  mission  requirements  given 
no  knowledge  of  the  time  of  hostilities  initiation,  or  location  of 
targets  at  this  time.  That  is,  the  mission  could  be  accomplished  with 
a  random  "starting  time"  for  hostilities,  and  with  the  placed  weapons 
also  in  random  positions  within  their  orbits  or  latitudes.  There  has 
been  no  attempt  to  relate  the  positions  of  the  targets  and  weapons  in 
time. 

Functional  Description 

Before  any  description  of  the  basic  HELBASE  functions  can  be  of 
value,  several  terms  used  frequently  must  be  defined. 

Basic  Definitions 

1.  Basing  Mode:  This  term  reflects  the  essence  of  HELBASE, 

and  refers  to  weapon  locations  on  the  ground,  in  aircraft,  or  in  space¬ 
craft.  The  mean  surface  of  the  earth  represents  a  single  basing  mode. 
The  airspace  above  the  mean  surface  of  the  earth  constitutes  the  second 
basing  mode.  The" third  and  last  mode  is  actually  a  continuous  set  of 
alternatives,  each  representing  a  user  selected  circular  orbit  altitude. 

2.  Weapon:  One  HEL  from  a  specific  weapon  type. 

3.  Weapon  Type:  A  specific  combination  of  a  set  of  HEL  weapon 
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characteristics  and  a  basing  mode.  A  set  of  weapon  characteristics 
Include  b sr.*  power,  beam  wavelength,  maximum  length  of  time  for  each 
firing  of  a  weapon,  minimum  time  between  firing  operations  of  a  weapon, 
maximum  number  of  times  each  weapon  is  operated,  maximum  firing  range, 
beam  waist  size,  and  a  probability  of  weapon  failure.  A  given  set  of 
the  weapon  characteristics  Is  combined  with  a  basing  mode  to  form 
a  weapon  type.  A  note  about  the  "maximum  firing  range"  is  apropos 
here.  This  quantity  is  not  a  limit  defined  by  the  inability  of  the 
weapon  to  negate  a  given  target  at  a  distance  past  this  range  as  a 
result  of  insufficient  power  applied  to  the  target.  Instead,  it  is  a 
rough  limit  offered  to  the  user  as  a  bound  on  the  weapon's  ability 
to  acquire.  Identify,  and  track  the  type  of  target  chosen  by  the  user 
for  that  weapon.  If  these  subsystems  do  not  constrain  the  use  of  the 
weapon,  the  default  value  may  be  accepted. 

4.  Weapon  System:  The  collection  of  weapons  selected  by  HELBASE 
from  the  available  weapons  types,  along  with  specific  locational  para¬ 
meters  for  each  weapon. 

5.  Target:  A  single  specific  satellite  upon  which  a  weapon  may 
be  trained  and  fired. 

6.  Target  Type:  A  group  of  targets  which  are  related  by  a 
common  set  of  mission  and  orbital  parameters. 

7.  Target  System:  The  collection  of  all  selected  target  types. 
For  example,  a  target  type  consisting  of  four  navigational  satellites 
could  form  a  target  system. 

8.  Mission:  The  operational  capability  around  which  the 
weapon  system  is  to  be  defined.  Target  type  selections  together  with 

mission  data  serve  to  fully  define  the  "mission." 
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9.  "Most  Efficient":  As  used  to  describe  the  functions  of 
HELBASE,  "most  efficient"  refers  to  the  weapon  system  (collection  of 
weapons  in  specific  orbits  or  on  specific  latitudes)  "built"  by  HELBASE. 

It  is  important  to  repeat  that  the  model  is  intended  to  arrive  only  at 
a  "first  cut"  approximation  of  the  "best"weapon  system.  Therefore,  the 
term  "most  efficient"  is  to  be  viewed  in  a  first  cut  approximation  sense. 
Without  doubt,  there  are  other  systems  of  weapons  that  are  more  efficient 
However,  their  definition  is  beyond  the  intent  of  this  first  cut  model. 

User  Inputs.  The  four  categories  of  user  inputs  are  Target  Types, 
Weapon  Types,  Mission  Data,  and  Battle  Management  Data.  Each  of  these 
areas  will  be  defined  in  terms  of  its  variables  and  the  effects  it  has 
on  the  functions  of  HELBASE.  Table  I  depicts  a  summary  of  all  inputs 
with  the  associated  default  values.  Most  of  the  information  in  this  table 
is  also  applicable  to  the  detailed  analysis  of  Chapter  III. 

As  defined  earlier,  a  target  type  is  a  group  of  targets  which  are 
related  by  a  common  set  of  mission  and  orbital  parameters.  For  example, 
two  communication  satellites  ir,  synchronous  orbit  may  constitute  a  target 
typo.  For  this  developmental  HELBASE  model,  the  user  is  restricted  to  the 
selection  of  any  combination  of  four  target  types  which  were  constructed 
to  be  reasonable  representations  of  the  range  of  "real  world"  possibilities. 
The  user  may  select  any  of  fifteen  (15)  different  combinations  of  target 
types  to  form  the  target  system.  There  is  currently  no  provision  for 
inputting  a  target  type  of  the  user's  design  (see  Figure  1). 

As  defined  above,  we  need  both  weapon  chararteri sties  and  a  basing 
mode  to  specify  a  weapon  type.  A  basing  mode  involves  one  selection 
(ground,  air  or  space)  and  (for  space)  one  additional  specification  for 


each  basing  mode  desired.  Selection  of  the  ground  mode  limits  the  deployment 
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TABLE  I 


Input  Category 
Target  Types 
Weapon  Types 


Mission  Data 


Battle  Management 


Dimensions  of  User  Inputs 


Variable  Default* 


0,  2.  3*  4)  * 

Wavelength  * 

Rowe*  * 

Wai  Size  1  M 

Total  Firing  Cycles 
Maximum  Firing  Time 
Recycle  Time  0 

Maximum  Range  «*> 

Probability  of  Failure  0 

Basing  Mode  * 

Altitude  0  Ground 

12  KM  Air 
*  Space 

Absorption  Coefficient  * 


Target  Types  * 

Target  Type  Negation 

Percentage  .  8 

Target  System  Negation 

Time  * 

Target  Type  Priority  1 

Data 

Target  Type  * 

Target  Type  Firing 

Priorities  4 

Irradiation  Time  Sorting  0 


Dimension 


Nanometers 

Watts 

Meters 


Seconds 

Seconds 

Kilometers 


Kilometers 


Hours 


♦Note:  An  asterisk  indicates  a  mandatory  entry 


TYPE  I!  SIX  SATELLITES  IN 
250  m  SLIGHTLY  ELLIPTICAL 
ORBITS 


TYPE  II!  TWO  SATELLITES  IN 
ASYNCHRONOUS  ORBITS 


TYPE  III!  THREE  SATELLITES 
IN  500  NM  CIRCULAR  ORBITS 


TYPE  IV!  THREE  SATELLITES 
IN  EQUATORIAL  SYNCHRONOUS 
ORBITS 


fig  1 .  Target  Types 
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of  weapons  with  the  associated  weapon  characteristics  to  any  location 
on  the  mean  surface  of  the  earth.  That  is,  HELBASE  will  deploy  the 
weapons  across  latitudes  so  as  to  meet  the  mission  requirements  with 
the  minimum  number  of  weapons  possible.  Selection  of  the  air  mode 
will  cause  deployment  of  weapons  in  aircraft  at  39,370  feet,  again 
across  latitudes.  Selection  of  an  aircraft  altitude  was  judged  to 
be  of  minimal  importance  regarding  laser  beam  attenuation,  and  so  the 
nominal  altitude  of  39,970  feet  (12  KM).  Selection  of  the  space 
basing  mode  requires  the  specification  of  a  circular  orbit  altitude 
in  addition  to  weapon  characteristics. 

As  a  demonstration  of  multiple  weapon  type  selections,  assume  a 
user  has  selected  the  ground  mode  with  weapon  characteristics  set  #1 
and  the  space  mode  with  weapon  characteristics  set  #2  along  with  an 
orbital  altitude  of  400  nautical  miles.  A  third  selection  is  the 
space  mode  with  weapon  characteristics  set  #2  (again)  along  with  an 
orbital  altitude  of  200  nautical  miles.  HELBASE  would  then  select 
the  "most  efficient"  combination  of  weapons  from  all  weapon  types  to 
satisfy  the  mission  requirements.  This  efficient  combination  may 
include  some  from  each  weapon  type,  or  aVI_  from  just  one  of  the  weapon 
types.  The  "most  efficient"  weapon  system  may  include  one  weapon  with 
weapon  characteristics  of  set  #1  at  55°  North  latitude,  no  weapons  at 
200  nautical  miles,  and  4  weapons  with  weapon  characteristics  from 
set  #2  at  400  nautical  miles. 

Mission  Data  inputs  define  the  mission  parameters  under  which 
the  selected  target  types  must  be  negated.  A  mission  is  not  simply  a 
matter  ''f  "melting  all  the  targets  in  the  shortest  time  possible"  but 
rather  "mission  denial"  (enough  radiation  on  enough  satellites  to 
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render  the  target  type  incapable  of  performing  its  mission)  in  a  speci¬ 
fied  maximum  length  of  time.  Toward  this  end,  the  user  inputs  Target 
Type  Negation  Percentages,  Target  System  Negation  Time,  and  Target  Type 
Priorities. 

The  Target  Type  Negation  Percentage  is  the  percentage  of  a  parti¬ 
cular  target  type  which  must  be  negated  to  insure  mission  denial.  For 
example,  we  may  have  to  negate  six  out  of  nine  satellites  of  a  particular 

target  type  (or  .67)  to  insure  tha>;  the  target  type  could  no  longer 

perform  its  mission.  A  target  type  negation  percentage  must  be  input 

for  each  target  type. 

The  Target  System  Negation  Time  is  the  maximum  time  the  user  will 
allow  the  weapon  system  to  accomplish  the  mission.  The  target  system 
negation  time  may  also  be  conceptualized  as  the  "maximum  mission  accom¬ 
plishment  time  allowed." 

Target  Type  Priorities  are  relative  weights  assigned  to  each  target 
type,  and  are  used  to  bias  the  deployment  of  the  weapon  system  toward 
those  target  types  with  the  higher  priorities.  That  is,  the  resulting 
weapon  system  locational  parameters  (ground-and  air-based  latitude,  space 
based  orbital  parameters)  and  to  a  lesser  degree  the  selection  of  weapon 
types  will  be  biased  in  favor  of  negating  the  higher  priority  target  types 
sooner.  The  priority  weights  are  values  of  1  to  1000.  Note  that  this  is 
a  strategic  concept,  in  that  the  deployment  of  the  weapons  is  affected. 
This  must  not  be  confused  with  target  type  firing  priorities,  which  is  a 
tactical  concept  discussed  below. 

Battle  Management  Data  inputs  reflect  the  tactical  concerns  of  the 
user.  Options  included  in  this  HELBASE  model  are  Target  Type  Firing 
Priorities  and  Irradiation  Time  Required.  (Other  options  not  included 
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are  listed  under  recommendations  in  Appendix  B).  The  user  may  select 
either,  neither,  or  both  options. 

If  the  user  wishes  to  use  the  engagement  tactic  of  assigning  target 
type  firing  priorities,  a  priority  of  one  is  assigned  to  the  highest 
priority  target  type  and  so  on  to  the  lowest  priority  target  type.  (An 
equal  priority  may  be  assigned  to  more  than  one  target  type).  Assign¬ 
ment  of  these  priorities  will  cause  relative  inefficiencies  in  the  use 
of  weapon  power,  but  will  tend  to  negate  higher  priority  target  types 
first.  That  is,  at  a  given  instant  in  time,  a  specific  weapon  may  have 
a  choice  from  many  tasget  opportunities.  The  target  type  firing  prior¬ 
ities  will  cause  this  "target  opportunity  set"  to  be  separated  into 
groups.  Further  separation  by  irradiation  time  required  is  possible 
within  each  of  these  groups,  if  desired  by  the  user. 

This  time  is  a  function  of  distance  to  the  target  and  target  hard¬ 
ness.  Obviously,  the  target  requiring  the  shortest  irradiation  time 
would  be  negated  first.  If  target  type  firing  priorities  are  not  also 
assigned^all  the  targets  in  the  target  opportunity  set  are  ordered  by 
least  irradiation  time.  If  priorities  are  also  assigned,  e  h  priority 
group  of  targets  is  ordered  by  least  irradiation  time. 

To  recap,  the  user  may  select:  1)  neither  option  (random  selection 
from  the  target  opportunity  set);  2)  either  option  (simple  ordering  by 
target  type  or  by  irradiation  time);  3)  £r  both  options  (target  type 
sorting  followed  by  irradiation  time  sorting  within  target  types). 

Functional  Flow 

•  A  broad  understanding  of  HELBASE  functional  flow  may  be  gained 
by  examination  of  figure  2.  Again,  the  simple  purpose  of  the  model 
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Fig  2.  HELBASE  Functional  Flow 
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is  to  translate  user  inputs  of  targets,  weapons,  tactics,  and  the  opera¬ 
tional  mission  into  a  first-order  approximation  of  the  "most  efficient" 
constellation  of  weapons  required  to  perform  the  mission.  "Most  efficient" 
can  be  restated  as  the  "minimum  number"  of  weapons  required  to  perform 
the  mission.  Again,  just  as  the  term  "most  efficient"  is  based  on  a 
first  order  approximation,  the  term  "minimum  number"  is  also  approximate. 

An  alternate  approach  to  understanding  of  HELBASE  functional  flow 
is  presented  in  figure  3  .  The  submodel  SPATIAL  is  just  the  collection 
of  SPHERE,  TARGET,  and  SIGHT,  all  of  which  determine  the  three-dimensional 
spatial  relationships  needed.  The  submodel  OPTIMIZE  acts  as  the  controll¬ 
er  of  HELBASE,  directing  all  model  activities  toward  the  definition  of 
the  "most  efficient"  weapon  system.  The  purposes  of  PROP  (for  propagrte) 
and  EFFECT  should  be  apparent.  BATTLE  allows  the  weapons  placed  by 
OPTIMIZE  to  fire  against  the  selected  target  system  over  the  mission 
negation  time. 

Basic  Physical  Principles 

The  fundamental  aspects  of  orbital  mechanics,  laser  propagation  and 
laser  effects,  as  used  in  HELBASE,  are  described  below. 

Orbital  Mechanics.  As  stated  in  Chapter  I  as  an  assumption,  no 
perturbations  in  orbital  behavior  are  considered.  The  earth  is  represent¬ 
ed  by  its  mass  located  at  the  center  of  a  sphere  of  radius  3443.923 
nautical  miles  (NM),  or  6378.145  kilometers  (KM).  Distances  are  referred 
to  in  either  NM  or  KM,  but  the  model  uses  KM  only.  In  the  same  light, 
earth  satellite  position  is  treated  as  the  classic  two  body  problem,  with 
perfect  elliptical  orbits. 

The  coordinate  reference  system  used  is  the  geocentric-equatorial 
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system,  which  has  its  origin  at  the  earth's  center  and  is  basically  fixed 
with  respect  to  the  stars.  This  choice  allows  us  to  easily  investigate 
the  behavior  of  ground-based  devices  whose  locations  need  only  be  repre¬ 
sented  by  latitudes.  Figure  4 depicts  this  reference  system.  Combining 
the  assumptions  above  (simple  two-body  mechanics)  with  this  coordinate 
system  provides  straightforward  relationships  between  the  target  system 
and  the  weapon  system.  Satellites  (both  target  and  weapon)  will  have 
"stationary"  orbits  relative  to  one  another,  and  weapons  displayed  on 
the  ground  or  in  aircraft  will  tend  to  rotate  along  latitudes  below. 

The  uses  of  this  reference  system  and  orbital  mechanics  are 
restricted  primarily  to  the  random  placement  of  satellites  into  orbits, 
and  the  movement  of  the  satellites  along  these  orbits  through  time. 

The  random  placement  of  a  satellite  into  an  orbit  means  that  five 
of  thu  six  orbital  elements  necessary  to  define  the  position  of  a 
sat'  'lite  are  held  constant  (see  figure  4  ),  and  the  sixth,  mean  anom¬ 
aly.  s  selected  randomly.  The  relationships  needed  for  this  random 
selection  are  as  follows  (Ref  5:185,  220-222). 

M  z  mean  anomaly.  M  varies  uniformly  from  0  to  2u. 
fl  T7.  eccentric  anomaly.  E  does  not  vary  uniformly,  but  as  a 
function  of  the  eccentricity  of  the  orbit. 

V— true  anomaly.  True  angle  from  perigee.  (See  figure  4  ). 

Pertinent  relationships  are: 

M  *  £  —  e  $  //7  £  cos  t  -  e  CO  j  H 

/  y-  c  cos  0 

Note  that  random  placement  of  satellites  into  their  orbits  is  equivalent 
to  randomizing  the  time  at  which  hostilities  begin. 
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co  argument  of  perigee 
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Fig  A.  Geocentric  Equatorial  Reference 


These  same  relationships  are  used  to  move  satellites  along  their 
respective  orbits  in  the  submodel  BATTLE. 

Laser  Propagation.  As  stated  in  the  introduction,  both  the 
propagation  and  laser  effects  submodels  form  the  lower  accuracy  bound 
on  the  performance  of  HELBASE.  (See  Chapter  V  for  methodology  suggested 
to  increase  the  effectiveness  of  HELBASE).  The  laser  propagation  sub¬ 
model  PROP  is  based  on  accountability  for  the  most  basic  two  of  the  many 
atmospheric  effects  on  a  laser  beam:  beam  expansion  (of  a  collimated 
beam),  and  atmospheric  attenuation  (for  the  ground  and  aircraft  modes 
only). 


For  propagation  of  a  beam  from  a  ground  -,  air  -,  or  space- 
based  HEL,  expansion  of  the  beam  will  occur  as  a  function  of  the  original 
"waist  size"  (Ref  7:17).  This  represents  the  initial  beam  radius 
(assumed  TEMoo  gaussian  beam  (Ref  8:94-97))of  a  non-focused  beam.  The 
expansion  of  the  beam  also  increases  with  distance  from  the  device  and 
with  wavelength  (Ref  7).  A  first  order  approximation  for  the  radius 
of  the  beam  of  wavelength\at  a  distance  2  from  the  weapon  is  given  by 

'4 

w  a>-  vJ  1+ 

l  VTfVlfcv  J  ,  where  2  is  "large" 


ard  w„  is  the  initial  beam  radius  at  the  waist  mentioned  above. 

Laser  beam  attenuation  (absorption)  by  the  atmosphere  can  be 
expressed  simply  as  a  form  of  Beer's  Law  (Ref  7).  We  ignore  the  added 
effects  of  particulates  and  aerosols,  accounting  only  for  the  decreasing 
density  of  the  atmosphere  as  the  beam  propagates  from  the  ground  or 
aircraft  toward  the  target. and  the  water  vapor  present  in  the  atmosphere. 
In  general ,  we  have 


,  -t' 

-X0e 


,  where  X  is  measured  along  a 
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direction  of  propagation  perpendicular  to  the  surface  of  the  earth,  N  (X) 

represents  the  molecular  density  (molecules/cm  )  of  the  atmosphere  as  a 

function  of  X,  andc'the  mean  molecular  cross  section  coefficient  for  the 

wavelength  considered  (Ref  8).  The  above  formula  may  be  approximated  by 
a  -jt/r.3* 


TN4  dx 

*X#  e  - 


This  approximation  is  not  accurate  for  particulates  and  aerosols,  but  it 
does  allow  the  calculation  of  beam  attenuation  without  finding  N  (X). 

Chapter  III  contains  a  continuation  of  the  derivation  of  the  forms 
the  above  formula  exhibits  for  ground  -  and  air  -  based  weapons.  Both 
attenuation  and  expansion  effects  are  combined  in  Chapter  III  in  order  to 
generalize  the  effect  on  beam  intensity  at  the  target. 

Laser  Effects.  Generally,  the  effect  (material  changes)  undergone 
by  a  target  when  irradiated  by  a  laser  beam  is  dependent  upon  the  proper¬ 
ties  of  the  material  being  irradiated,  the  beam,  and  the  physical  relation¬ 
ship  between  the  two. 

Material  properties  at  the  initiation  of  laser  irradiation  are 
absorptivity,  thermal  conductivity,  temperature,  specific  heat,  density 
and  latent  heats  of  fusion  and  vaporization  (Ref  7;  9:1-2).  Beam  proper¬ 
ties  of  interest  are  wavelength,  intensity  and  length  of  time  the  beam  is 
applied.  We  will  not  consider,  allow,  nor  model  "pulse"  lasers  of  such 
short  duration  and  high  power  so  as  to  cause  mechanical  pressure  on  the 
target.  Conduction  of  the  absorbed  heat  away  from  the  target  area  is 
important  in  the  case  of  longer  irradiation  times  and  lower  power  levels. 
Other  complexities,  such  as  the  thickness  of  the  material  and  the  effects 
of  melt  retention  and  vaporization  of  the  outer  surface  of  the  irradiated 
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area,  may  also  be  considered  for  a  specific  application  of  Interest. 

Methodology  for  determining  whether  a  HELBASE  target  Is  destroyed 
is  contained  in  Appendix  A. 

Submodel  Descriptions 

Each  submodel  in  HELBASE  is  functionally  analyzed  below.  Basic 
Interrelationships  with  other  models  are  described  along  with  the  sub¬ 
model  inputs  and  outputs.  In  some  cases  a  submodel  may  be  used  in 
slightly  different  ways  by  other  submodels.  The  three  submodels  of 
SPATIAL,  (that  is,  SPHERE,  TARGET,  and  SIGHT),  will  be  followed  by  the 
two  submodels  which  deal  with  the  laser  beam,  PROP  and  EFFECT.  SIDEN 
(for  sighting  density)  which  is  a  part  of  the  next  model,  OPTIMIZE, 
and  BATTLE  finish  the  list. 

SPHERE.  The  purpose  of  SPHERE  is  tn  calculate  and  store  the 
cartesian  coordinates  of  1650  points  which  define  a  "sphere  of  operations" 
(See  figure  5).  These  1650  points  are  spaced  approximately  5°  apart  in 
latitude  and  longitude,  and  thus  approximate  a  5°  uniform  distribution 
of  points  over  the  sphere.  This  sphere  of  operations  represents  one  of 
the  basing  modes  input  by  the  user.  A  400  nautical  mile  (740.8  kilometers) 
circular  orbit  basing  mode  is  represented  by  a  sphere  of  operations  of 
radius  3844  NM  (7119  KM).  The  1650  points  are  then  spread  uniformly 
over  this  sphere  of  operations,  and  each  point  represents  an  instant¬ 
aneous  weapon  location  possibility.  Great  circles  around  this  sphere 
therefore  represent  possible  weapon  orbits. 

It  is  important  to  state  here  that  this  sphere  of  operations 
is  fixed  in  the  geocentric  equatorial  coordinate  system,  so  that,  over 
time,  the  earth  turns  beneath  the  sphere.  This  restricts  alternatives 
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Fig  5*  Submodel  SPHERE 
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In  the  ground-based  mode  to  latitude  choicest  and  the  spacing  of  our 
1650  points  In  fact  restricts  our  accuracy  to  5°  latitudinal  Increments. 

So  SPHERE  calculates  and  stores  the  cartesian  coordinates  of 
all  1650  points  on  a  sphere  of  operations,  which  may  represent  the 
operational  geometry  of  ground-,  air-,  or  space-based  weapons.  A  sphere 
of  operations  must  be  defined  for  each  basing  mode  selected.  SPHERE 
requires  as  Inputs  the  radius  of  the  sphere  of  operations,  and  outputs 
the  cartesian  coordinates. 

TARGET.  While  SPHERE  generates  the  operational  geometry  of  weapon 
locational  possibilities,  TARGET  does  something  a  little  different  for 
the  target  satellites.  TARGET  simply  places  the  targets  into  their 
proper  orbits  in  a  random  manner,  and  stores  the  cartesian  coordinates 
of  all  target  locations.  The  TARGET  algorithm  involves  a  random  select¬ 
ion  of  a  mean  anomaly  from  a  uniform  (0,  2t()  distribution,  a  root  solut¬ 
ion  to  Kepler's  equation  for  the  eccentric  anomaly,  and  then  a  simple 
conversion  to  the  true  anomaly,  (See  fig  6).  Since  TARGET  stores  the 
orbital  elements  for  all  four  allowed  target  types,  it  only  requires  an 
input  of  which  target  types  were  selected.  The  output  is  the  cartesian 
coordinates  of  each  target  in  the  target  system  (selected  target  types). 
TARGET,  SPHERE,  and  two  submodels  not  yet  addressed  are  used  by  the  next 
submodel,  SIGHT. 

SIGHT.  Before  we  can  discuss  the  purpose  of  SIGHT,  we  need  to 
define  a  "sighting"  as  used  here.  Rather  than  the  common  understanding 
of  an  open  line-of-sight  between  a  point  and  the  target,  this  "sighting" 
carries  with  it  information  about  the  feasibility  of  negating  the  target 
from  the  point  in  question.  If  we  indeed  have  a  feasible  opportunity  to 
negate  the  target,  then  we  score  a  "1".  If  any.  of  the  variables  which 
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Ffg  6.  Submodel  TARGET  Logic  Diagram 
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must  be  checked  to  insure  a  feasible  "sighting"  is  a  value  which  prohibits 
the  negation  opportunity,  then  we  score  a  "0".  Items  checked  to  determine 
whether  we  have  a  "sighting"  (feasible  negation  opportunity)  are  inter¬ 
ference  of  the  earth,  maximum  irradiation  time  restriction  (for  weapon 

power  conservation),  and  maximum  effective  range.  (See  definition  of  range 
in  "basic  definitions"  section). 

To  add  a  little  more  information  to  our  binary  "0,  1"  score  we 
divide  by  the  irradiation  time  required  to  negate  the  target  if  the  score 
is  1.  This  will  allow  a  more  efficient  selection  of  weapon  locations  by 
OPTIMIZE.  See  Chapter  II  for  the  development  of  this  attribute.  So 
SIGHT  provides  weighted  information  concerning  the  feasibil ity  of  a 
negation  opportunity  for  given  target  and  weapons  location.  A  value  of 
0  indicates  that  no  opportunity  exists,  and  increasing  values  indicate 
increasing  power  efficiency  of  feasible  shots.  See  figure  7  for  the 
SIGHT  logic  diagram.  Required  information  is  the  location  and  type  of 
the  target,  and  location  and  type  of  the  weapon.  SIGHT  will  then  provide 
the  target  sighting  efficiency  value  as  described  above.  SIGHT  is  used 
by  SIDEN  (a  part  of  OPTIMIZE)  and  BATTLE. 

PROP.  This  submodel,  along  with  EFFECT,  forms  the  lower  bound  on 
HELBASE  accuracy.  Depending  on  the  use  of  HELBASE  intended  by  the  user, 
both  propagation  and  laser  effects  modeled  in  PROP  and  EFFECT  may  need 
to  be  validated  by  the  user  and  possibly  expanded.  (See  Chapter  V  for 
details).  However,  for  many  uses  the  propagation  and  effects  models 
described  here  are  sufficient.  (See  Chapter  III  for  formulation  and 
supporting  analyses). 
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Fig  7.  Submodel  SIGHT  Logic  Diagram 
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PROP  simply  takes  the  wavelengths,  initial  beam  intensity,  beam 


waist  size,  and  the  geometry  between  the  weapon  and  the  target,  and 
translates  this  information  into  an  intensity  of  the  beam  at  the  target. 
Algorithms  are  included  for  ground  to  space,  air  to  space  and  space  to 
space  beam  spreading  and  attenuation. 

EFFECT.  In  order  to  determine  the  irradiation  time  required  to 
negate  a  target  f0r  which  the  beam  intensity  has  been  calculated  by 
PROP,  we  use  the  submodel  EFFECT.  This  model  relies  on  a  data  base 
discussed  in  Appendix  A.  Assumptions  and  limitations  of  transforming 
this  data  to  be  used  with  an  incoming  wavelength,  are  discussed  in 
Appendix  A.  With  this  data  base,  EFFECT  uses  the  beam  intensity  at 
the  target  (and  its  wavelength)  to  calculate  the  irradiation  time 
required  to  negate  the  target. 

SIDEN.  The  purpose  of  SI  DEN  is  to  generate  an  "average  target 
system  sighting  efficiency"  value  for  each  point  on  each  sphere  of 
operations.  It  does  this  essentially  by  using  TARGET  to  randomly  place 
the  targets  into  their  orbits  and  then  using  SIGHT  to  calculate  the 
target  system  sighting  efficiency  of  each  point  on  each  sphere  of  oper¬ 
ations.  This  implies  that  SIGHT  must  be  called  for  each  target  for  one 
given  point,  and  this  process  then  repeated  for  all  the  points  in  that 
sphere  and  then  for  all  spheres  of  operations.  This  information  is 
retained  and  the  target  system  drawn  is  thrown  out  and  a  new  one  drawn. 
Target  system  sighting  efficiencies  are  again  calculated  and  again  new 
target  system  locations  are  drawn.  This  process  is  repeated  until  an 
"average  target  system  sighting  efficieny"  value  is  reached  for  each 
point  on  each  sphere  of  operations.  The  average  target  system  sighting 
efficiency  value  then  represents  the  percentage  of  all  targets  which  can 
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be  "seen"  from  that  point  (on  the  average),  weighted  by  the  inverse  of 
the  irradiation  time  required.  Chapter  III  provides  the  supporting 
analysis  for  this  construct,  and  explains  the  rationale  to  use  this 
value  as  a  basis  for  weapon  placement  by  OPTIMIZE. 

The  following  explanations  may  help  to  clarify  the  meaning  of  this 
important  efficiency  attribute. 

One  weapon  location  looking  at  one  target:  Negation 
opportunity  value  or  target  sighting  efficiency. 

One  weapon  location  looking  at  all  targets  in  the  system 
(as  stochastically  placed  by  TARGET):  Target  system 
sighting  efficiency. 

One  weapon  location  looking  at  all  targets  randomly  drawn 
many  times:  Average  target  system  sighting  efficiency. 

These  average  target  system  sighting  efficiencies  are  used  for 
the  deployment  decision  (by  OPTIMIZE)  of  only  one  weapon.  Once  this 
weapon  is  placed,  the  efficiency  values  are  no  longer  used.  If  another 
weapon  is  needed  to  satisfy  mission  requirements,  a  new  efficiency  must 
be  calculated  for  each  point.  The  new  efficiency  values  incorporate  the 
effect  of  the  weapon  (s)  already  placed  by  using  as  the  target  system 
the  remaining  targets  at  the  end  of  a  BATTLE  (as  opposed  to  using  a 
complete  target  system  as  generated  by  TARGET)  in  the  calculation  of  the 
target  sighting  efficiencies.  Figure  8  is  a  logic  diagram  reflecting  the 
functional  description  above.  SIDEN  is  a  major  part  of  the  OPTIMIZE 
submodel . 

OPTIMIZE.  As  mentioned  above,  OPTIMIZE  is  the  controller  of 

HELBASE,  and  has  as  its  primary  function  the  selection  of  the  "most 

efficient"  placement  of  weapons  in  the  basing  modes  selected  by  the  user. 

It  does  so  by  the  simple  aggregation  of  the  efficiency  attribute  "average 
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Submodel  SI  DEN  Logic  Diagram 


target  system  sighting  efficiencies"  (discussed  in  the  last  section) 
around  orbits  (or  latitudes  for  the  case  of  the  ground  or  air-based 
spheres  of  operations)*  and  then  selection  of  the  orbit  or  latitude 
with  the  maximum  aggregated  value  of  average  target  system  sighting 
efficiencies.  Details  of  this  selection  procedure  are  provided  in 
Chapter  III. 

In  addition  to  placing  the  weapons,  OPTIMIZE  checks  whether  the 
weapon  system  consisting  of  all  placed  weapons  can  meet  the  mission 
requirements.  This  check  is  made  by  running  a  BATTLE  for  the  mission 
negation  time  specified  in  mission  data.  If  the  mission  requirements 
are  met  by  the  end  of  the  mission  negation  time,  the  weapon  system  is 
considered  sufficient  to  meet  the  user's  needs.  This  check  is  made 
many  times  in  order  to  develop  some  confidence  in  the  ability  of  the 
weapon  system  to  accomplish  the  mission.  (Confidence  interval  statis¬ 
tics  in  Chapter  IV)  (See  figure  9). 

BATTLE.  The  only  submodel  which  advances  over  time  is  BATTLE, 
which  uses  the  weapon  system  supplied  by  OPTIMIZE  to  engage  the  target 
system  over  the  mission  negation  time.  In  this  submodel  tactics 
characterized  by  battle  management  data  are  used. 

The  targets  and  weapons  are  randomly  placed  in  their  orbits  or 

latitudes,  the  BATTLE  clock  is  set  to  zero,  and  the  engagement  is  begun 

within  the  constraints  imposed  by  battle  management  data.  The  tactics 

available  to  the  user  are  "target  type  firing  priorities"  and'irradiation 

time  required."  The  user  options  are  fourfold:  1)  Let  both  tactics 

default  (this  would  result  in  random  target  selection)  2)  Assignment 

of  target  type  firing  priorities  would  cause  targets  to  be  ordered  in 

terms  of  priority  (target  selection  within  priorities  would  be  random) 

3)  The  tactic  of  selecting  targets  by  least  to  most  irradiation  time 
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Fig  9-  Submodel  OPTIMIZE  Logic  Diagram 


required  without  priorities  would  seem  to  be  the  most  energy  conservative 
of  the  four  available  tactics,  4)  Last,  the  user  may  elect  to  sort  the 
targets  by  priority  and  order  by  least  Irradiation  time  within  each 
priority. 

The  above  tactics,  as  selected  by  the  user,  are  applied  at  an 
instantaneous  time  during  che  progress  of  BATTLE  (See  figure  10).  This 
procedure  incorporates  a  certain  degree  of  approximation  concerning 
the  irradiation  time  versus  the  unallowed  movement  of  both  target  and 
weapon  during  this  time.  The  effect  of  this  procedure  is  minimized 
by  keeping  the  BATTLE  clock  advancement  time  on  the  order  of  some  average 
irradiation  time,  say  about  five  minutes.  See  Chapter  III  for  a  detailed 
analysis  of  BATTLE  clock  advancement  time  calculation.  At  t  -  0  on  the 
BATTLE  clock,  BATTLE  must  call  SIGHT  for  all  targets  over  all  weapons. 

The  irradiation  times  required  tor  those  shots  that  are  feasible  are 

/ 

recorded.  The  selected  tactics  are  then  employed  to  let  the  weapons 
fire  against  the  targets.  The  BATTLE  clock  is  then  advanced,  and  all 
remaining  targets  and  weapons  are  advanced  a  corresponding  distance 
along  their  orbits  (or  along  their  latitudes).  This  process  is  repeated 
until  the  BATTLE  clock  reaches  the  mission  negation  time.  OPTIMIZE  then 
examines  the  results  of  the  BATTLE. 

This  chapter  is  concluded  with  a  table  summarizing  the  HELBASE 
submodels,  and  their  functions,  inputs  and  outputs.  (See  table  II). 
Chapter  III  forms  the  logical  and  mathematical  rationale  for  the  HELBASE 
functions  described  above. 


START 


fig  10.  Submodel  BATTLE  Logic  Diagram  (1  of  2) 
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Submodel  BATTLE  Logic  Diagram  (2  of  2) 
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TABLE  II 


Summary  of  Submodel  Functions 


SUBMODEL 

FUNCTION 

INPUTS 

OUTPUTS 

SPHERE 

Calculate  and  store  cartesian 
coordinates  of  1550  points  on 
a  sphere  of  operations. 

Radius 

Cartesian  coordinates 
>f  all  points. 

TARGET 

Place  targets  randomly  into 
orbits. 

Target  type 
selections. 

Cartesian  coordinates 
>f  all  targets. 

SIGHT 

Calculates  target  sighting 
efficiency. 

Target  and 
weapon  loca¬ 
tion  weapon 
type  and  tar¬ 
get  type. 

Target  sighting  effi 
iency. 

PROP 

Calculates  beam  intensity 
at  the  target. 

Weapon  type, 
weapon  and 
target 
locations. 

Beam  intensity  at 
target. 

EFFECT 

Calculates  irradiation  time 
required  to  negate  a  target. 

Target  type 
and  beam 
intensity  at 
target. 

Irradiation  time 
required. 

SI  DEN 

Establishes  the  target  system 
sighting  efficiency  for  each 
point  on  each  sphere  of  opera¬ 
tions. 

Weapon  type, 
target  type, 
spheres  of 
operation. 

Values  of  target 
system  sighting 
efficiency. 

OPTIMIZE 

Select  weapon  orbits  or  lati¬ 
tudes  and  check  mission 
requirements. 

Average  tar¬ 
get  system 
sighting 
efficiencies. 

Weapon  system 
locational  para¬ 
meters. 

BATTLE 

Engage  targets  over  mission 
negation  time  in  accordance 
with  battle  management  data. 

Weapon  system 
target  system 
mission  data, 
battle  manage 
ment  data. 

Number  of  targets 
negated. 
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Ill  MODEL  DEVELOPMENT 


This  chapter  contains  the  background  analyses  upon  which  various 
parU  vi  the  model  are  based.  These  derivations  of  various  concepts  and 
methodologies  are  Included  to  support  the  functions  of  the  model.  Infre¬ 
quently,  methodological  attempts  which  failed  may  be  discussed  as  back¬ 
ground  material  for  those  developments  which  were  useable. 

Sphere  of  Operations 

The  concept  of  a  "sphere  of  operations,"  which  is  a  stationary 
sphere  in  the  celestial  reference  system,  is  used  to  renresent  all  possible 
locations  of  a  particular  weapon  type  (see  figure  11  ).  It  is  a  good  approx¬ 
imation  for  ground-based  weapons,  slightly  less  good  for  aircraft-based 
weapons,  and  the  most  restrictive  for  space-based  weapons. 

Ground  Basing  Mode.  For  ground-based  weapons,  the  sphere  of 
operations  is  the  surface  of  the  earth.  A  ground-based  weapon  describes 
a  latitudinal  circle  around  its  sphere  of  operations  as  the  earth  rotates 
360°  around  its  axis.  Therefore,  this  path  becomes  a  decision  alternative 
(a  possible  laser  weapon  location)  for  the  submodel  OPTIMIZE  (fig.  9  ). 

The  continuous  set  of  latitude  choices  is  restricted  to  the  discrete  set 
of  latitudes  spaced  at  5°  intervals  (see  figure  12)  for  simulation  pur¬ 
poses. 

Aircraft  Basing  Mode.  A  similar  idea  is  used  for  aircraft- 
based  weapons,  although  the  concept  of  latitudinal  paths  is  not  as  accur¬ 
ate,  The  aircraft  is  restricted  to  "orbiting"  a  ground  reference  point 
during  weapon  operation,  so  that  its  path  also  approximates  a  latitudinal 
circle  as  the  earth  turns  beneath  the  sphere  of  operations.  Obviously 
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12.  Ground-Based  Alternative  Latitudes 
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an  aircraft  traveling  at  450  knots  could  cause  some  deviation  from  the 
"latitudinal  circular  path"  concept,  particularly  at  higher  latitudes. 

This  restriction  on  aircraft  movement  is  conservative  in  relation 
to  the  weapon  system  needed  to  meet  the  mission  requirements.  Chapter  V 
discusses  the  analytical  improvement  of  weapon  system  performance  which 
may  be  obtained  by  removing  this  restriction  after  the  simulation  has 
placed  the  full  weapon  system. 

Space  Basing  Mode.  A  sphere  of  operations  for  a  space-based 
weapon  is  defined  by  a  given  radius  from  the  earth's  center.  A  continuous 
set  of  circular  orbit  possibilities  is  approximated  by  825  orbital  possi¬ 
bilities  spaced  approximately  uniformly  over  the  surface  of  the  sphere 
of  operations.  Each  of  these  825  orbital  alternatives  is  defined  by  a 
unit  momentum  vector,  each  pointing  from  the  center  of  the  sphere  of 
operations  to  one  of  the  825  points,  (fig  13  ). 

Placement  of  Points.  The  original  intent  was  tc  refine  a  uniform 
spherical  distribution  of  points,  each  a  distance  of  five  arc  degrees 
from  its  nearest  neighbors.  After  a  search  for  and  an  attempted  develop¬ 
ment  of  an  exact  three  dimensional  function  revealed  the  limitations  of 
the  class  of  regular  polyhedrons,  an  approximation  was  calculated  as 
follows. 

For  each  latitude,  beginning  with  the  equator  as  the  first  latitude 
and  proceeding  North  in  five  degree  increments,  a  "nearest  latitude"  was 
found  which  allowed  an  integer  number  of  five  degree  arcs  (whose  length 
was  defined  at  the  equator)  to  be  spaced  around  its  circumference.  For 
example,  the  equator,  at  0°  latitude,  has  a  circumference  of  360°,  allow¬ 
ing  it  to  be  divided  into  exactly  72  5°  arcs.  However,  the  15°  North 
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latitude  circumference  is  only  347.7333°  in  circumference,  when  measured 
in  terms  of  the  circumferential  distance  a  1  arc  degree  would  describe 
at  the  equator  (great  circle).  Therefore,  only  69.5467  5°  arcs  could  be 
transcribed  around  its  circumference.  However,  if  we  move  South  to  the 
13.5362°  latitude,  exactly  70  5°  arcs  can  be  defined  around  the  circum¬ 
ference. 

This  technique  of  finding  the  "nearest  latitude"  was  repeated  over 
all  latitudes  from  5°  North  to  the  northern  pole  of  the  sphere.  The 
results  are  shown  in  table  III  .  The  southern  hemisphere  is  a  mirror 
image  of  the  northern  hemisphere.  Figure  11  depicts  a  typical  sphere  of 
operations. 


TABLE  HI 

Number  of  "Great  Circle  Five  Degree  Arcs" 
for  North  Latitudes  of  a  Sphere  of  Operations 


True  Latitude 

Nearest  Latitude 

Number  of  5° 

0° 

0° 

72 

5° 

5° 

71.72602 

10° 

9.5603 

71 

15° 

13.5362 

70 

o 

o 

CsJ 

19.1881 

68 

25° 

25.4744 

65 

30° 

30.5584 

62 

35° 

34.9708 

59 

o 

o 

*3- 

40.1918 

55 

45° 

44.9005 

51 

50° 

50.2910 

46 

46 


True  Latitude 

Nearest  Latitude 

Number  of  5°  Arcs 

55° 

55.2885 

41 

60° 

60.0000 

36 

65° 

65.3757 

30 

70° 

69.6825 

25 

75° 

74.6991 

19 

80° 

80.4059 

12 

00 

tn 

o 

85.2198 

6 

o 

o 

CTi 

90.0000 

(1  PT) 

Submodel  SPHERE  calculates  the  cartesian  coordinates  of  all  of 
these  points  and  stores  them  for  use  later  in  the  model  (fig  14).  See 
Appendix  C  for  a  FORTRAN  listing  of  this  submodel. 

Submodel  TARGET 

TARGET  uses  the  orbital  parameters  of  the  satellite  targets 
selected  to  randomly  place  them  into  their  orbits.  The  only  "phase 
relationships  allowed  in  this  first  HELBASE  version  is  between  the  sub- 
synchronous  and  synchronous  targets. 

The  algorithm  used  is  the  "mean  anomaly  to  eccentric  anomaly  to 
true  anomaly  technique"  (Ref  5:182-185).  A  mean  anomaly  is  randomly 
drawn  from  a  unifjrm  (0,  2'iT)  distribution.  The  mean  anomaly  M  is  re¬ 
lated  to  the  eccentric  anomaly  E  by: 

jvl-  t-  e.s./>£ 

where  e  is  the  eccentricity  of  the  orbit.  A  root  solution  numerical 
technique  is  used  to  find  the  zero  value  (EQ)  of: 
f  (E)  -  E  -  esinE  -  M, 

where  e  and  M  are  known.  The  technique  used  (ZBRENT,  International 
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Mathematical  and  Statistical  Libraries),  requires  that  E  values  (initial 
guesses)  be  input  such  that  their  functional  values  (F  (E))have  differ¬ 
ent  algebraic  signs.  The  simple  algorithm  of  making  the  first  guess  E^  = 
M  -  2e  and  the  second  guess  =  M  +  2e  guarantees  that  f(E^)<0  and 
f(E2>0. 

The  EQ  value  is  then  used  to  calculate  the  true  anomaly  by: 

-i)  =  COS~‘(  ^ 

\  e  cos  tc-l  / 

Note  that  the  root  solution  procedure  is  unnecessary  when  e  is 
zero,  since  in  this  case  the  mean,  eccentric,  and  true  anomalies  are  the 
same. 

Sighting  Ef fi ci ency  and  Submodel  SIGHT 

"Sighting  Efficiency"  is  not  only  the  basis  of  submodel  SIGHT, 

it  is  the  fundamental  attribute  upon  which  submodel  OPTIMIZE  selects 

the  "most  efficient"  weapons  to  construct  the  final  weapon  system. 

Several  attributes  were  considered  before  the  decision  was  made  to 

employ  "sighting  efficiency."  Each  of  the  alternatives  had  drawbacks, 

* 

but  the  attribute  selected  seems  to  encompass  the  essence  of  "most 
efficient"  and  has  validity  for  a  majority  of  possible  model  users. 

Development  of  Sighting  Efficiency.  The  objective  under  consider¬ 
ation  here  is  "to  place  the  next  weapon  in  the  most  efficient  position." 
Efficiency  is  defined  in  terms  of  determining  the  weapon  system  composed 
of  the  least  number  of  weapons  needed  to  accomplish  the  mission.  A 
definition  refinement  process  was  used  in  determining  the  attribute  which 
could  in  a  reasonable  way  measure  the  attainment  of  this  objective. 

The  efficiency  attribute  was  to  be  defined  first  for  a  point  (1  of 
the  1650)  on  the  sphere  of  operations,  and  then  extended  to  measure  the 
efficiency  of  the  weapon  path  (latitudinal  circle  or  orbit). 
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The  first  attribute  considered  was  simply  the  average  percent- 
age  of  the  total  target  system  which  could  be  "sighted"  at  any  point  in 
time.  The  approach  was  to  randomly  distribute  the  satellites  in  the 
target  system,  and  then  to  determine  whether  each  of  them  could  be 
"seen"  from  each  point  on  the  sphere  of  operations.  A  "sighting"  in 
this  case  was  defined  as  1)  satisfaction  of  non-interference  of  the 
earth  with  the  shot,  2)  a  check  of  the  maximum  weapon  range  (a  function 
of  pointing  accuracy  and  3)  a  check  of  the  maximum  irradiation  time 
allowed  against  what  was  required  to  negate  the  target.  A  binary  cod¬ 
ing  system  was  to  be  used  to  record  the  satisfaction  of  these  checks. 
Non-interference  with  the  earth,  less  than  maximum  range  and  less  than 
maximum  irradiation  time  would  be  coded  as  a  "1",  and  violation  of  any 
of  these  three  constraints  would  be  coded  as  a  "0".  Adding  and  then 
averaging  these  values  over  many  random  draws  of  the  targets  would 
yield  the  desired  average  percentage.  For  instance,  the  sighting  of 
5,  6,  12  and  2  targets  out  of  15  over  4  stochastic  draws  of  the  targets 
would  yield  a  value  of  .4167  for  that  point  on  a  specific  sphere  of 
operations.  That  is,  an  average  of  about  42%  of  the  targets  could  be 
seen  from  that  point  at  a  given  time. 

The  impact  of  the  weapons  already  placed  into  the  weapon  system 
by  the  submodel  OPTIMIZE  was  represented  by  discounting  any  "signtirg" 
which  could  also  be  "sighted"  by  any  of  the  weapons  already  placed. 

This  was  intended  to  bias  tne  placement  .of  the  next  weapon  away  from 
the  sphere  of  operations  areas  already  "covered"  by  placed  weapons. 

This  attribute  was  soon  recognized  to  be  too  simple  and  incomplete. 

It  was  recognized  that  this  measurement  paid  no  attention  to  the 
irradiation  time  required  to  negate  the  target  (other  than  an  upper  limit 
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check).  This  meant  that  a  negation  which  would  take  15  seconds  from  a 
range  of  10,000  KM  and  a  negation  of  2  seconds  at  1000  KM  would  both 
receive  a  value  of  "1".  Therefore,  the  submodel  OPTIMIZE  would  as  likely 
pick  a  point  characterized  by  longer  irradiation  times  as  one  character¬ 
ized  by  shorter  times.  So  this  first  attribute  attempt  violated  the 
"completeness"  property,  ( Ref  1 0: 50)  in  that  it  proved  inadequate  to  meas¬ 
ure  the  degree  to  which  the  objective  was  obtained.  It  was  apparent  that 
"most  efficient"  points,  and  then  weapon  paths,  should  be  measured  by 
something  more  than  a  binary  (0,  1)  "no  sighting"  or  "sighting"  criteria. 

The  second  attempt  at  the  definition  of  an  attribute  by  which  to 

measure  the  "most  efficient"  point  on  the  sphere  of  operations  therefore 

included  the  irradiation  time  required  for  a  specific  weapon  placed  at  that 

point  to  negate  a  target.  This  was  accomplished  by  dividing  the  previous 

binary  "1"  by  the  "irradiation  time  required"  to  negate  the  target.  For 

example,  a  sighting  of  2  out  of  10  targets  at  negation  times  of  10  seconds 

apiece  would  yield  a  value  of  1_  +  1  =  .20*  Also,  the  sighting  of  1  tar- 

10  10 

get  out  of  10  at  a  negation  time  of  5  seconds  would  yield  1  =  .20.  These 

5 

events  are  therefore  considered  to  be  of  equal  value.  The  definition  of  a 
"sighting"  remained  the  same  as  before:  line-of-sight  interference  with  the 
earth,  maximum  range,  and  maximum  irradiation  time  allowed. 

At  this  point  it  was  noted  that  discounting  al 1  sightings  of  targets 
which  could  also  be  seen  by  any  of  the  weapons  already  placed  tended  to  over¬ 
compensate  for  the  weapons  placed.  No  accounting  for  the  capability  of  the 
placed  weapon  to  actually  negate  the  targets  it  could  "sec"  was  being  made. 

It  was  feasible  that  a  placed  weapon  with  only  one  operating  cycle  in  its 
design  lifetime  could  prohibit  the  allocation  of  any  other  weapons  to  its 
portion  of  the  sphere  of  operations. 
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Allowing  one  discounting  of  a  sighting  by  a  placed  weapon  seemed 
too  conservative,  so  the  manner  by  which  the  effect  of  the  weapons 
already  placed  could  be  represented  became  a  problem.  The  matter  was 
finally  resolved  by  taking  as  the  initial  target  system  (to  be  looked 
at  by  all  points  over  all  spheres  of  operations)  the  targets  which  re¬ 
mained  after  a  BATTLE. 

The  submodel  BATTLE  was  to  be  run  for  the  input  mission  negation 
time.  Any  targets  left  after  this  time  would  be  looked  at,  and  the 
"sightings"  scored  by  one  divided  by  the  "irradiation  time  required" 
as  before.  The  placement  of  the  next  weapon  would  therefore  be  made 
contingent  upon  the  effect  of  all  weapons  previously  placed. 

An  additional  benefit  occurred  when  it  was  decided  to  take  the 
remaining  target  set  after  reaching  the  user  input  mission  negation 
time  in  a  BATTLE.  Additional  information,  in  the  forms  of  total  weapon 
operating  cycles  (or  total  irradiation  time  allowed),  and  weapon  re¬ 
cycle  time  were  implicitly  included  in  the  effect  the  "already  placed" 
weapons  would  have  on  the  target  system. 

Submodel  SIGHT.  The  submodel  SIGHT  was  designed  specifically  to 
provide  a  value  of  sighting  efficiency,  for  a  "one-on-one"  look  for  a 
point  on  a  sphere  of  operations  at  a  target,  (see  fig  7  ).  If  a  line- 
of-sight  check  did  not  prohibit  a  "shot"  (see  figures  15,16,&l7),then 
checks  would  be  made  for  exceeding  the  maximum  allowable  irradiation 
time  for  a  weapon  operation,  and  the  maximum  distance  allowable  (function 
of  pointing  accuracy). 

Calculation  examples  are  shown  in  the  section  which  discusses 
the  submodel  "SIDEN."  See  Appendix  C  for  a  listing  of  SIGHT  in  Fortran. 


CASE  I!  GOOD  LINE  OF  SIGHT  (PROJI*Re<,  EOT  COM  <  ^*/2) 


TARGET 


COM 

PROJI  IS  THE  DISTANCE 
FROM  THE  EARTH  S  CENTER 
PERPENDICULAR  TO  THE 
TARGET/WEAPON  POSITION 
VECTOR. 


COM  IS  THE  ANGLE  BETWEEN 
THE  TARGET/WEAPON  VECTOR 
AND  THE  VECTOR  FROM  THE 
EARTH  S  CENTER  THROUGH  THE 
WEAPON. 


CASE  II !  BAD  LINE  OF  SIGHT  (PROJI<Re<>  COM  /2) 


PROJI 


PROJI 


Fig  15.  Line  of  Sight,  Ground-Based  Case 
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CASE  l:  GOOD  LINE  OF  SIGHT  (PROJl>Re  ,  COM<fi/2) 


Submodel  SIDEN.  The  function  of  SIDEN  is  to  collect  the  data 


necessary  and  then  to  calculate  the  "target  sighting  efficiency"  for  each 
point  in  each  sphere  of  operations. 

For  one  point  on  a  sphere  of  operations  looking  at  one  target,  we  have 


k  * 


S  *  1  ,  where 
t 


S  =  Sighting  efficiency,  and 
t  =  Irradiation  time  required  to  negate  the  target. 

If  we  extend  this  concept  to  n  targets  (in  the  target  system),  where 
1,  2,  ...»  n,  then  _  A  A 

••Is*  “SX 

K*l  KM 


where  t^  is  the  irradiation  time  required  to  negate  the  k  th  target,  and 
s  is  the  average  of  S k  over  all  n  targets. 

Extending  the  "sighting  efficiency"  concept  to  multiple  stochastic 
draws  of  the  targets  (j  *  1,  2,  ...jn)  for  m  iterations,  we  have 

m  n 


EFF 


jM  KM 


where 


EFF  =  Average  sighting  efficiency  of  the  point  over 
m  stochastic  draws  of  n  targets. 

=  Negation  time  required  for  the  jth  drawing  of 
the  kth  target. 

This  average  sighting  efficiency  for  a  point  on  the  sphere  of  operations 
is  extended  to  include  the  path  average  sighting  efficiency  in  the  submodel 
OPTIMIZE  (See  Appendix  D), 

Submodel  PROP 

This  section  investigates  the  basis  and  development  of  the  submodel 
PROP  (propagation)  whose  purpose  is  to  translate  beam  intensity  at  the  source 
to  beam  intensity  at  the  target. 
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The  only  effects  modeled  in  this  first-order  model  (see  Chapter  II, 
Limitations)  are  beam  spreading  (for  all  basing  modes),  and  beam  absorption 
(for  air-or  ground-based  weapons).  Beam  spreading  does  not  include  thermal 
blooming,  and  beam  absorption  does  not  include  aerosols  or  particulates. 

Beam  Absorption.  Absorption  of  the  laser  beam  during  its  transition 
through  the  atmosphere  is  modeled  on  a  standard  atmosphere.  The  funda¬ 
mental  formulation  is  an  extension  of  Beer's  law  to  a  non-homoceneous 
medium. 

For  an  aircraft-based  weapon  firing  vertically  up  through  the 
atmosphere,  the  intensity  at  the  target  is  given  by  ( Ref 1 7) : 


•^74  ~  e 


fj 


-x/S.S' 

N*T  e  dx 

h-y 


where 


Ht/*=  Beam  intensity  at  the  target, 
le-  beam  intensity  at  the  weapon, 
h =  Weapon  height  from  the  surface  of  the  earth, 
A/  =  Molecular  density  of  the  air, 

<T  =  Absorption  coefficient, 

=  Vertical  height  from  weapon  to  target,  and 
X  =  Distance  al;  .-i  beam  from  weapon  to  target. 


(See  figure  18 ). 

If  the  effect  of  shooting  at  an  angle  (  6  )  from  the  local 


vertical  is  included  in  our  beam  absorption,  we  have 

or 

j7A=i„ec«eJ"-< 


d> 


Note  that  for  a  ground-based  weapon,  h  -  y  =  0,  so  that 
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(Ground-Based  Weapon) 

A  simple  approximation  can  be  made  to  avoid  wasting  computer 
time  on  numerical  integrations  (PROP  may  be  called  on  the  order  of 
one  million  times  during  a  run  of  HELBASE).  A  simple  sensitivity 
analysis  showed  that  very  little  difference  in  the  value  of  It  (less 
than  .0001%)  is  found  if  the  upper  limit  of  integration  is  set  at  200  KM. 

2CG 


Therefore,  we  have 


osG  Jh-y 


NT 

CO- 


or  approximately 


—  *TeJ  & 


-(W-Y) 
-6.S  NTCi 
CO3  0 


Beam  Spreading.  An  independent  phenomenon  associated  with  laser 
beam  transmission  over  long  distances  is  beam  divergence,  or  beam  spread¬ 
ing.  The  fundamental  formulation  of  the  beam  divergence  is  based  upon  a 
collimated  (non-focused)  beam  of  waist  size  wQ  at  the  source.  A  collimated 
beam  will  spread  proportionately  to  its  wavelength  and  inversely  proport¬ 
ionately  to  its  waist  size  (Ref  7;  6:17).  The  "waist  size"  (a  gaussian 
beam  is  assumed)  can  be  thought  of  as  the  radius  from  beam  center  to  a 
point  at  which  the  irradiance  is  1/e  of  that  at  beam  center  (Ref  6  :247). 

The  waist  size  of  the  spreading  beam  can  ("far"  from  the  source)  be  express¬ 
ed  as: 
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where 


o>c  *  Waist  size  at  the  source, 


tc'g  -  Waist  size  at  target, 

A  =  Beam  wavelength,  and 
2  =  Distance  from  weapon  to  target  (Ref  7). 

Simple  geometry  provides  the  relationship  for  the  reduction  in  beam 
intensity  resulting  from  beam  spreading 


(  3-t 5  =  Intensity  at  target  reduced  by  beam  spreading) 


lTS=Io(f7)J 


Combined  Beam  Spreading  and  Absorption.  Combining  the  results 
of  the  last  two  sections  can  be  expressed  as: 


where  IQ  has  been  replaced  by  1^ 


the  attenuated  beam  intensity  at  the  target.  This  relationship  expands 


to 


e 


-&.S  N' 
an  e 


-  (h-y)  /  8- S 


(-0 


This  expression  of  the  final  beam  intensity  at  the  target  is  appropo 
for  ground  and  air-based  weapons.  For  space-based  weapons,  the  simple 
beam  expansion  formula  is  used. 

Submodel  EFFECT. 

The  submodel  EFFECT  is  designed  to  "use"  data  of  the  general  type 
found  in  the  Foreign  Technology  Division/Lockheed  Missile  and  Space 
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Company’s  "Project  ESS"  analyses.  For  a  look  at  representative  data  from 
these  studies,  see  Appendix  A. 

The  data  can  be  reduced  (Appendix  A)  to  a  quadratic  expression  of 
the  form 

at2  +  bi  +  (c  -  JjT  jt  )  - o 

where 


a,  b,  c,  d,  are  quadratic  curve  fitting  and  data  constants  based  on 


the  data, 

t  =  The  negation  time  (irradiation  time)  required  to  negate  the 
target), 

=  Wavelength  of  the  beam  under  consideration,  and 
\  =  Beam  intensity  at  the  target. 

This  relation  can  be  solved  for  the  negation  time,  which  is  the  output 


of  EFFECT: 


i=-b+  [  It)] 


'4 


Zo. 

This  simple  algorithm  may  be  inspected  in  Appendix  C. 


Submodel  OPTIMIZE 

The  most  important  function  of  the  submodel  OPTIMIZE  is  to  identify 
the  single  weapon  path  (latitudinal  circle  or  orbit)  which  is  the  "most 
efficient,"  because  onto  this  path  the  next  weapon  allocated  will  be 
placed. 

Various  techniques  were  investigated,  some  of  which  may  be  the 
basis  of  the  future  for  model  expansion  (see  Appendix  B).  The  methodology 
finally  selected  involves  a  straightforward  check  of  all  the  possible  paths, 
with  selection  of  the  maximum  as  the  "most  efficient." 


iStSiM  UHj  1 
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For  each  ground  -  or  air-based  weapon,  a  "path  sighting  efficiency" 
around  each  latitude  across  the  sphere  of  operations  (in  5°  intervals)  is 
calculated.  This  amounts  to  37  possible  paths  for  each  weapon.  For  each 
space-based  weapon,  there  are  825  possible  orbits.  These  orbits  are  de¬ 
fined  by  the  direction  of  the  angular  momentum  vectors  (see  fig  11)  going 
through  the  upper  hemisphere  of  points  defining  a  sphere  of  operations. 

Calculating  a  "path  sighting  efficiency"  for  each  of  these  825 
orbits  meant  finding  the  locus  of  points  on  the  sphere  of  operations 
which  most  closely  approximated  the  orbital  path.  However,  once  the  se¬ 
quence  of  points  constituting  each  orbital  path  were  identified,  they 
could  be  used  for  all  space-based  spheres  of  operation.  The  idea  was  to 
find  the  36  points  which  would  define  each  of  the  825  orbits,  and  to  store 
this  data  until  needed  (permanent  storage).  For  this  purpose,  a  sub¬ 
routine  called  ORB  was  developed  (see  Appendix  C). 

Subroutine  ORB.  The  subroutine  ORB  was  needed  only  to  define  the 
36  points  on  each  of  the  825  orbits.  A  fundamental  part  of  this  subroutine 
is  the  capability  to  rotate  the  three-dimensional  cartesian  coordinate 
axes.  A  brief  development  of  these  axes  rotations  is  presented  here  prior 
to  looking  at  the  methodology  of  determining  the  36  points  of  each  orbit. 

Figure  19  shows  a  typical  vector  in  both  the  spherical  and  cartesian 
coordinate  systems. 

Transformation  of  spherical  coordinates  to  cartesian  coordinates  is 
a  one-to-one  transformation,  and  is  simply 

X=  psinG  cor3  q 
p  Sifl0  5irt 
2r  pCCS  0 
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Transformation  of  cartesian  coordinates  to  spherical  coordinates 
is  not  a  one-to-one  transformation.  We  may  say  that 

rn — 5 - —  ^  .  -i  / 

?=VX  +H  ^  Sin  l 

/  q  \  »  and 

w(tJ 

where  the  following  algebraic  sign 

s 

combination  exist: 

jt  -fc 

(*-, <4~ /Zt  )  ( *1 ‘i-  z ' )  ( *+) (a t ft, z -  )  (*+.  T.Z-) 

For  the  rest  of  the  octants,  we  may 
say  that  cos”1  (a/(>s*r>0) 

where,  the  following 
algebraic  sign  combinations  exist: 

(>w z+  X*  •  <  r)  (*■;  ht.  )0c  h".  ) 

♦Note:  For  these  algebraic  sign 

combinations,  ,  * 

=  -s in'  [h /%*»'&) 

To  rotate  the  cartesian  coordinate 


reference  system  through  an  angle  0  about  the  Z  axis,  we  take  (Ref  A:  77-80) 


To  rotate  this  reference  system  through  an  angle  (  Q-Ti/Z  )  about 


the  y  axis,  we  take 

Cos  (Q-Ti/z) 
Hz  ~  O 

&.  Sin  CQ  -  n  /z) 


-sin  (e-  h/2)1  fXi' 
o  H, 

cos  (  6-  Ti/z)  2> 


Therefore  to  rotate  the  reference  system  through  an  angle  0  about 
the  Z  axis,  and  then  through  an  angle  (0  -u/2)about  the  rotated  y  axis, 
we  can  say  that 


rccs(<3-fi/2l  0  -  sin  (o-ti/z\  "lrco^ 

Hi.  -  O  *  O  -sincfi 

2x.  si*lG-fi/£l  O  cosCe-fi/tl  o 


Sm<JS  0~lp£ 

-Sincf>  C0$9  0 

_  O  O  i  J 1 2b 


Xil  "oc5  (  Q  -H/S^cos^  CCI  lo-fi/SK.'si'ifl  -S<n(G -lifer 

Wi|  =  -Stn<$  CC5<£  ^  4t> 


—  Stn<£ 

si<i  te-fi/sKccsqn 


Sin  CG'ii/2^S*n(f»  COiftH'/^UuSk 


Similarly,  a  rotation  matrix  was  developed  to  rotate  a  given  reference 
system  backward  through  (<Vz-0)  about  the  y  axis,  and  then  through  (-<£) 
about  the  rotated  2  axis 


Ho  - 


cos<}i  cos Cfi/2  -oN  -st/up 

Sin<|  cox  til/2  -  Gl  cos^ 

si»(n/z-Q)  o 


-CCi(p  Cm  £  ri/2-  Ci)  -<  -ti- 
-  sm  Sm  £  ff/«r- (31 

Ce£  Cit/2.-  6l  J  l?c 
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We  are  now  ready  to  investigate  the  methodology  used  to  determine 
the  points  used  in  each  of  the  orbits.  The  methodology  is  that  employ¬ 
ed  in  subroutine  ORB  (see  Appendix  C).  Notation  used  in  this  section 
is  as  follows: 

Cartesian  coordinates:  * -  Point  number 

* -  Reference  system 


Spherical  coordinates: 


Point  number 
Reference  system 


Axes: 


Xo 

Hr 


Original  x  axis 
Rotated  Z  axis. 


We  start  with  the  cartesian  and  spherical  coordinates  of  a  point 
on  the  sphere  of  operations  which  represents  the  unit  angular  momentum 
vector.  The  superscript  M  is  used  on  the  coordinates  to  designate  a  unit 
momentum  vector  as  opposed  to  a  point  number  on  an  orbit  (see  fig. 20  ). 

The  first  step  is  to  rotate  the  original  (Xo,  Yo,  Zo  reference 
system  through  the  angle  M  about  the  ZQ  axis,  and  then  through  the 
angle  (G^'i'/c  )  about  the  YR  axis  via  the  following  matrix: 


g" 

Zc 


cc~  CGc -rt/Slcc-f  ( ($>£*)  cos  (Oc  -  £./>  (Qc’ 


—  *tn 


CCS  (  ) 


O 


[_  cci  C  1  s*n  (  Gj'-i 1/t  \  «•*»  cci  (Q?~ 


Coordinates  of  the  momentum 
vector  in  the  rotated  reference 
frame. 


Xc 

Ho 

CJ 

'S 


Coordinates  of  the  momentum 
vector  in  the  original  reference 
frame. 
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ORB  Axis  Rotations 


These  rotations  place  the  XR  axis  along  the  momentum  vector,  and 
the  2^  axis  through  the  first  point  of  the  orbit. 

The  coordinates  of  this  first  point  on  the  orbit  are  ©<*-rt/2. 

The  task  now  is  to  express  the  coordinates  of  this  point  in  terms 
of  the  original  reference  system  (Xo^  ,  ,  ZQ^  ).  We  therefore 

must  rotate  the  coordinates  of  the  first  point  backward  through 
-  ( u/2.  )  about  the  YD  axis,  and  then  through  (-0  M)  about  the 

K  0 


Z^  axis, 
o 


If 

Hi 

»«■ 

m 

co: 

s.*»  (<$")  cus(\i/z-  Gj?\ 


6c  \ 


ces 
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-  cos ( <f0rt >  j.</» ( t7 /z-  ef\ 
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Hi 
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Similarly,  to  obtain  the  coordinates  of  the  jth  point  of  the  ith 
orbit,  add  5°  to  the  0  of  the  (j  -  1)  st  point,  and  solve  for  the  cartes¬ 
ian  coordinates  in  the  rotated  reference  frame.  We  know  that 

^  s  ^  s  constant, 

=r  it fz-  constant,  and 
Gg  -  0^+5° 

.  * 

Xg  —  {> cci  (  siX  (  i-0 

=  fc C4  (6c> 

We  can  then  solve  for  the  (  of 

the  original  reference  system. 

The  last  step  in  the  process  of  defining  the  points  which  describes 


a  given  orbit  over  the  surface  of  a  sphere  of  operations  is  to  find  the 
closest  point  to  our  above  precisely  defined  orbital  point.  The  algorithm 


for  accomplishing  this  is  shown  in  Appendix  C. 

Submodel  BATTLE 

The  purpose  of  the  submodel  BATTLE  is  to  allow  the  placed  weapons 
to  fire  against  the  targets  over  the  mission  negation  time,  under  the 
constraints  imposed  by  the  user.  FiglOin  Chapter  II  shows  the  function¬ 
al  flow  of  BATTLE,  and  the  Fortran  subroutine  is  reproduced  in  Appendix 
C. 

The  functional  description  in  Chapter  II,  Table  I  ,  and  Appendix 
C  describe  the  various  combinations  of  tactics  available  to  the  user. 

This  section  concentrates  on  two  modeling  aspects  of  this  many-on-many 
engagement  submodel. 

The  first  matter  to  be  investigated  is  the  selection  of  the  "proper" 
t  for  use  in  advancing  the  BATTLE  clock.  Items  to  be  considered  in  the 
selection  of  a  At  are:  1)  the  angular  travel  of  the  satellite  which  has 
the  highest  angular  velocity  at  perigee,  2)  the  mission  negation  time, 

3)  the  minimum  of  all  weapon  recovery  times,  and  4)  computational  cost. 

BATTLE  At.  The  amount  of  an  orbit  covered  by  the  targets  and  the 
orbital  weapons  during  a  At  should  heuristically  be  kept  small  enough  to 
allow  target  acquisition  by  weapons  which  could  "sight"  the  target  during 
its  perigee  passage.  The  probable  target  types  with  the  highest  eccentric- 
ity  and  the  longest  semi-major  axis  will  have  these  approximate  character¬ 
istics:  e  =  .90  and 
a  9400  KM. 

From  the  expression  for  the  incremental  time  necessary  to  traverse 
an  incremental  portion  of  the  true  anomaly  (Ref  5:31), 

di^CrVM  do 
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where 


r  =  Instantaneous  radius, 
h  *  Angular  momentum,  and 
=  True  anomaly, 
we  can  say  (for  "small'&t  and/W),  that 

/It  s 

F  • 

From  h=  [&(  i-e1-) 

where  5.  cJf6  x  I0y  KM*/52^,  ^  0>-  a(i-el 


we  can  say  that  (for  satellites  near 


perigee),  «. <  i - <s\  J/kti  0/  (i-e^ 

So  for  the  case  of  e  =  .91, 
a  =  9356.96  KM,  and  At>  =  15°  (.2618  radians. 

At  *  7.33  seconds. 


So,  we  are  considering  a  lower  bound  for  BATTLE  At  of  about  7  seconds. 
The  second  item  to  be  considered  is  the  user  defined  mission  negation  time. 

A  short  mission  negation  time,  on  the  order  of  minutes  as  opposed 
to  hours,  may  create  irregularities  in  functions  of  BATTLE  if  the  BATTLE 
At  clock  advancement  increment  becomes  a  "large"  portion  of  the  mission 
negation  time.  For  example,  a  mission  negation  time  of  1  hour,  a  BATTLE 
At  of  5  minutes,  and  an  average  target  irradition  time  of  10  seconds  would 
allow  only  20  opportunities  for  the  weapons  to  fire  against  the  targets. 

A  shorter  At  would  allow  more  firing  opportunities,  and  a  closer  approxi¬ 
mation  of  realizing  the  full  potential  of  the  weapon  system.  So  we  must 
decrease  the  BATTLE  At  or  institute  some  lower  limit  on  the  user  input 


mission  negation  time.  For  example,  mission  negation  time  divided  by  At 
should  not  be  too  much  less  than  the  maximum  firing  time  allowed. 
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Weapon  recycle  tine  is  a  factor  in  that  a  BATTLE  At  which  is  much 
larger  than  the  shortest  recovery  time  would  ^llow  the  creation  of  in¬ 
efficiencies  in  weapon  employment  which  are  purely  the  result  of  the 
modeling  process. 

A  pervasive  item  which  is  a  re*l  limit  on  all  simulations  is  the 
computer  time  (cost)  required  to  accomplish  the  simulation.  A  BATTLE 
At  which  is  shorter  than  that  required  by  the  above  three  items  wastes 
computer  processing  time.  This  constraint  then  would  have  us  push  the 
BATTLE  At  to  as  large  a  value  as  possible. 

Reviewing  the  constraints  on  the  BATTLE  At,  we  have  lower  bounds 
formed  by:  1)  the  angular  travel  of  the  "fastest"  satellite  during  the 
At,  2)  a  fraction  of  the  mission  negation  time,  and  3)  the  least  recovery 
time,  and  an  upward  pressure  caused  by  computational  cost  consideration. 

Mission  negation  times  from  an  hour  to  several  days  are  consider¬ 
ed  practical.  The  7  second  lower  bound  based  on  the  angular  distance 
traveled  by  the  target  satellite  with  the  greatest  eccentricity  and  semi- 
major  axis  can  be  substantially  ignored  with  only  a  minor  impact  on  BATTLE 
performance.  Not  only  is  the  satellite  near  its  perigee  a  very  small 
portion  of  the  orbital  period,  it  is  also  a  more  demanding  pointing 
and  tracking  problem  to  actack  a  target  going  through  perigee  for  most 
weapon  configurations.  Letting  the  BATTLE  At  become  much  greater  than 
weapon  recycle  time  will  tend  to  cause  some  inefficiencies  in  weapon 
employment  when  weapons  which  are  useable  (have  been  recycled)  cannot  be 
used  because  the  BATTLE  At  has  not  been  reached.  However,  the  impact  of 
this  inefficiency  is  reduced  by  keeping  the  translation  of  the  satellites 
along  their  orbital  paths  down  to  a  reasonable  size,  so  that  the  targets 
will  generally  still  be  available  at  the  end  of  the  BATTLE  At. 
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A  At  of  5  minutes  would  allow  a  At?  of  about  .3  radians  or  17  degrees 
for  a  typical  target  from  the  target  system  allowed.  A  BATTLE  At  of  much 
more  than  five  minutes  would  introduce  discontinuities  in  the  progression 
of  the  many-on-many  engagements.  Unfortunately,  5  minute  At's  over  a 
2  day  mission  negation  time  would  mean  576  iterations  of  the  BATTLE 
functions. 

However,  a  BATTLE  At  of  five  minutes  is  the  original  value  used. 

A  sensitivity  analysis  of  model  output  to  the  value  of  the  BATTLE  /It  is 
certainly  warranted  (see  Chapter  IV). 

Progression  of  the  Satellites  Around  Thei r  Orbits  at  Each  At. 

Both  the  target  satellites  (through  subroutine  TARGET)  and  the  weapon 
satellites  are  assigned  an  initial  mean  anomaly.  The  effect  on  the  mean 
anomaly  of  a  time  increment  At  is  given  by 

AM  a  r.  At  (Ref  A :  1 85 )  where  n  =^T.  For  each  satellite,  then, 
a  new  mean  anomaly  M-,  =  Mq  +AM  can  be  calculated,  and  the  root  solution 
procedures  employed  in  TARGET  can  be  employed  to  calculate  the  positions 
of  all  the  weanons  and  satellites  at  t  +  At. 

This  chapter  has  investigated  the  underlying  physical  processes, 
methodologies,  and  techniques  for  the  complicated  submodels  of  HELBASE. 
The  detailed  information  herein  must  be  viewed  in  the  context  of  Chapter 
I  and  II,  which  ..;afine  the  constraints  on  and  functional  flow  of  each  of 
these  submodels. 
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IV.  Verification  and  Validation 


The  basic  methodology  proposed  by  Haley  and  Ghelber  (Ref  A)  was 
employed  in  the  HELBASE  verification/val idation  process.  This  effort 
(A  Methodology  for  Validation  of  Complex  Multi -Variable  Military  Computer 
ized  Models)  is  well  based  in  contemporary  thought  concerning  the  veri¬ 
fication  and  validation  of  military  simulation  models,  and  draws  on  the 
work  of  Finge  / Naylor,  Schlesinger,  Tytula,  and  others  (Refl2,13,14). 

In  this  approach,  model  verification  (assurance  that  the  model  functions 
as  intended)  is  a  part  of  the  overall  validation  process. 

The  process  is  divided  into  four  parts:  1)  conceptual  validity  -- 
an  assessment  of  the  purpose,  desired  accuracy,  assumptions,  limitations 
and  model  structure,  2)  verification  --  using  several  contemporary 
techniques,  3)  credibility  —  a  combination  of  face  validation  (expert 
judgment)  and  sensitivity  analysis,  and  4)  confidence  —  ways  to  help 
instill  some  degree  of  confidence  about  the  model  in  the  user  or  decision 
maker. 

Conceptual  Val idity 

Assessment  by  a  potential  user  of  the  conceptual  validity  of  HELBASE 
must  be  derived  from  Chapters  I  and  II.  The  model  purpose,  degree  of 
accuracy,  assumptions,  limitations  and  functional  structure  are  described 
in  detail  in  those  chapters. 

Verification 

This  part  of  the  overall  validation  process  is  primarily  concerned 
with  the  "mechanical"  validity  of  HELBASE,  and  is  composed  of  four  parts: 
1)  a  structured  walk-through,  2)  verification  of  technical  physical 
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processes,  3)  simulation  of  predictable  states,  and  4)  testing  of 
stochastic  events. 

Structured  Walk-Through.  This  part  of  verification  establishes 
the  event-path  integrity  of  the  model,  and  was  accomplished  for  the 
applicable  submodels  (modules)  as  well  as  for  the  overall  HELBASE  model. 

A  structured  walk-through  to  check  the  proper  treatment  of  all 
legal  variable  combinations,  was  accomplished  for  the  following:  INPUT, 
TARGET,  SIDEN,  OPTIMIZE,  and  BATTLE.  (See  Chapters  II  and  III  for 
functional  descriptions  and  the  mathematical  basis  of  each). 

The  treatment  of  the  INPUT  section  was  primarily  concerned  with  check¬ 
ing  the  default  mechanisms  for  weapon  types,  mission  data,  and  battle 
management  data.  In  all  cases,  the  user  must  input  a  zero  whenever  the 
model  default  value  is  wanted.  The  following  defaults  were  checked: 

Weapon  Types  —  Beam  waist  size, 

maximum  firing  cycles  allowed, 

maximum  firing  time  per  operating  cycle, 

recycle  time, 

maximum  range  allowed, 

probability  of  failure,  and 

weapon  attitude. 

Mission  Data  —  Target  type  negation  percentage,  and  target  type 
priority. 

Battle  Management  Data  -  Firing  priorities,  and  sorting  by  irradiation 
time. 

Submodel  TARGET  has  three  basic  event  paths  which  were  checked.  The 
normal  path  for  the  random  selection  of  target  locations  is  from  a  randomly 
selected  mean  anomaly  to  the  corresponding  eccentric  anomaly  and  finally 


to  the  true  anomaly.  This  path  was  checked  by  inspection  of  the  output 
and  by  hand  calculation  of  satellite  position  based  on  the  random  number 
drawn.  There  are  two  exceptions  to  this  standard  event  path.  First, 
target  types  two  and  three  (subsynchronous  and  synchronous  satellites) 
have  mean  anomalies  which  are  related.  The  first  satellite  of  the  two 
sub-synchronous  satellites  is  placed  randomly,  and  the  mean  anomaly  of  the 
second  satellite  is  180°  from  that  of  the  first  satellite.  Similarly, 
the  mean  anomalies  of  the  second  and  third  synchronous  satellites  are 
120°  on  either  side  of  the  randomly  selected  mean  anomaly  of  the  first 
satellite.  The  second  exception  to  the  normal  TARGET  event  Path  occurs 
when  the  target  orbital  eccentricity  is  zero,  as  for  target  types  three 
and  four.  In  these  cases,  the  mean  anomaly,  eccentric  anomaly  and  true 
anomaly  are  equal . 

We  have  now  investigated  the  event  path  integrity  of  submodels  INPUT 
and  TARGET.  SIDEN,  OPTIMIZE,  and  BATTLE  remain. 

Submodel  SIDEN  calls  submodel  SIGHT  for  "all"  targets  over  all  of  the 
points  on  each  sphere  of  operations.  The  italics  around  "all"  indicate 
the  dual  event  path  nature  of  this  submodel.  If  SIDEN  is  called  from  sub¬ 
model  OPTIMIZE,  the  target  system  selected  by  the  user  (with  all  selected 
targets  present)  is  looked  at  ("sighted")  by  each  o^  e  1650  points  on 
each  of  the  defined  spheres  of  operation.  However,  a  call  of  SIDEN  from 
the  submodel  BATTLE  causes  the  deletion  from  sighting  consideration  of  any 
targets  which  did  not  survive  the  BATTLE.  In  addition,  the  positions  of 
the  targets  are  randomly  drawn  by  TARGET  for  a  SIDEN  call  from  OPTIMIZE 
and  are  taken  as  they  are  at  the  end  of  a  BATTLE  for  a  call  from  BATTLE. 

The  two  basic  event  paths  in  the  "next  weapon  placement"  section  of 


75 


submodel  OPTIMIZE  are  concerned  with  searching  over  all  feasible  lati¬ 
tudes  for  ground-and  air-basing  modes,  and  over  all  feasible  orbits  for 
the  space-basing  mode.  The  two  event  paths  become  one  again  when  the 
path  matrix  is  searched  for  a  maximum  (see  Appendix  C). 

Event  path  integrity  for  the  submodel  BATTLE  is  a  little  more  complex. 
The  following  list  summarizes  the  decision  points  within  BATTLE: 

Initial  target  and  weapon  (random)  positioning, 
or  movement  along  paths  (latitudes  or  orbits), 

Calling  of  SIGHT  to  input  values  to  the  "battle 
efficiency"  (BATTLE  F)  matrix  -  call  SIGHT  if  the 
target  has  survived  thus  far  and  the  weapon  has  not 
met  its  total  operating  cycle  limit  or  is  currently 
firing  or  being  recycled. 

Firing  the  weapons  against  the  targets  —  the  user 
may  have  selected  (in  battle  management  data)  either 
firing  priorities  "or"  irradiation  time  sorting'1,  or 
both,  or  neither. 

Stopping  rule  option  -  BATTLE  may  be  stopped  by  the 
clock  reaching  the  mission  negation  time,  or  by  the 
negation  of  all  of  the  targets  prior  to  the  mission 
negation  time. 

Each  of  these  event  paths  was  checked  several  times.  An  aid  in  the 
verification  process  was  a  liberal  use  of  explanatory  comment  cards  before 
critical  decision  points  (see  Appendix  C). 

This  completes  the  structured  walk  through  of  the  HELBASE  sub¬ 
models.  The  second  part  of  the  verification  process  consists  of  a 
brief  review  of  the  HELBASE  processes  which  are  based  on  the  laws  of 
physics,  and  is  followed  by  a  look  at  predictable  model  output  and  a 
discussion  of  tests  applicable  to  HELBASE  stochastic  processes. 

Verification  of  Technical  Physical  Processes.  The  only  submodels 
which  reflect  complex  physical  processes  are  SPHERE,  TARGET,  PROP,  EFFECT, 
and  SIGHT. 

The  phenomena  underlying  each  of  these  submodels  was  covered  in 
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detail  in  Chapter  III.  In  addition,  hand  calculations  were  compared  to 
submodel  outputs  to  check  the  proper  coding  of  these  relationships. 

Simulation  of  Predictable  States.  The  HELBASE  model  purpose  of  placing 
weapons  in  a  "most  efficient"  manner  in  response  to  definition  of  weapon 
types,  a  target  system,  a  mission,  and  tactics  brings  to  mind  several 
combinations  of  inputs  for  which  an  output  weapon  system  could  be  predicted. 

Looking  at  the  weapon  type  input  variables,  several  outcome  states  can 

be  predicted  for  simple  combinations  of  inputs.  Definition  of  two  weapon 

types  in  a  space  basing  mode  in  an  exactly  identical  manner  with  the 

exception  of  wavelength  (and  the  corresponding  absorption  coefficient) 

should  lead  to  a  weapon  system  composed  entirely  of  weapons  of  the  type 

with  the  shorter  wavelength.  In  a  similar  way,  the  following  predictable 

states  should  occur  when  all  inputs  but  the  one  indicated  are  held  constant: 

Increasing  weapon  power  level  should  decrease  the  total 
number  of  weapons  required. 

Increasing  the  maximum  number  of  firing  cycles  allowed 
should  decrease  the  number  of  weapons  required,  unless  the 
mission  negation  time  is  of  the  same  order  of  magnitude  as 
the  number  of  operating  cycles  times  the  maximum  firing 
time  per  operating  cycle. 

Reducing  the  recycle  time  required  should  lead  to  fewer 
weapons  needed. 

Increasing  the  maximum  allowable  firing  range  should  lead 
to  fewer  weapons  unless  the  range  was  already  enough  to 
reach  all  or  most  of  the  targets. 

Combinations  of  certain  target  types  with  specific  weapon  types 
can  be  used  to  get  predictable  results.  Selection  of  type  four  (synchronous 
orbit)  targets  with  no  other  types  should  lead  to  placement  of  weapons  at 
or  near  the  equatorial  latitudes  and  equatorial  orbits.  Definition  of 
weapon  types  in  the  ground-  and  air-basing  modes  along  with  selection 
of  the  single  target  type  two  (subsynchronous)  should  lead  to  weapon 

placements  along  higher  North  latitudes. 
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Decreasing  the  target  type  negation  percentage  or  increasing  the 
target  system  negation  time  should  tend  to  decrease  the  number  of  weapons 
required  to  meet  the  mission  requirements.  Setting  one  target  type 
priority  to  anything  other  than  the  default  value  of  one,  while  leaving 
other  target  type  priorities  at  the  default  value  should  cause  inefficien¬ 
cies  in  weapon  placement  and  therefore  cause  an  increase  in  numbers  of 
weapons  needed. 

In  battle  management  data  inputs,  sorting  by  irradiation  time  should 
lead  to  the  used  fewer  weapons  than  the  default  random  target  selection. 
Similarly,  the  assignment  of  firing  priorities  as  opposed  to  acceptance 
of  the  default  value  should  lead  to  the  need  for  more  weapons  to  meet 
the  mission  requirements.  No  "predictable  state"  checks  have  been  made. 

Testing  of  Stochastic  Events.  The  only  random  selections  in  the  HEL- 
BASE  model  occur  in  submodels  TARGET  and  BATTLE,  and  are  from  the  simple 
uniform  (0,  2T»)  distribution.  The  simple  nature  and  use  of  these  random 
selection  procedures  does  not  seem  to  warrant  goodness-of-fit  tests  on  the 
random  selections  made. 

However,  there  are  three  areas  in  the  model  where  stochastic  analysis 
seem  appropriate.  The  first  concerns  the  "average  target  system  sighting 
efficiency  values  stored  in  the  SUMEFF  matrix  by  OPTIMIZE.  The  second 
area  is  the  derived  variable  "SUCCESS"  which  reflects  the  outcome  of  a 
BATTLE.  The  last  area  is  the  output  of  the  HELBASE  model  itself. 

As  a  brief  review,  the  values  stored  in  the  SUMEFF  (5  by  1650)  matrix 
are  the"average  target  system  sighting  efficiency"  values  corresponding  to 
the"efficiency"  of  a  single  point  on  a  specific  sphere  of  operations.  For 


example,  SUMEFF  (1,  1)  contains  the  "average  target  system  sighting 
efficiency"  value  for  point  1  on  the  sphere  of  operations  of  weapon  type 
#1.  Similarly,  the  value  stored  in  SUMEFF  (4,1650)  is  the  efficiency  of 
point  1650  on  the  sphere  of  operations  of  weapon  type  #4. 

Since  these  SUMEFF  values  are  used  by  OPTIMIZE  to  find  the  "most 
efficient"  path  on  which  to  place  the  next  weapon,  it  seems  reasonable 
to  want  these  values  to  be  good  estimates  of  their  "true"  values.  A 
"true"  value  in  the  SUMEFF  matrix  would  be  found  only  if  SIDEN  were  called 
infinite  number  of  times  by  OPTIMIZE.  Since  this  is  not  feasible,  esti¬ 
mates  of  the  "true"  values  in  SUMEFF  must  be  found  within  some  confidence 
interval  with  some  degree  of  confidence  that  the  true  value  actually  lies 
within  that  interval. 

Each  call  of  SIDEN  by  OPTIMIZE  produce  a  "target  system  sighting 
efficiency"  value  for  all  points  on  all  spheres  of  operation.  These 
values,  to  be  input  to  the  SUMEFF  matrix,  reflect  averages  of  the  efficiency 
of  a  given  point  overall  targets  (see  Chapter  III).  These  values  are 
averaged  into  the  SUMEFF  values  by  OPTIMIZE  (again,  see  Chapter  III  and 
also  Appendix  C).  This  is  why  a  SUMEFF  value  is  referred  to  as  an  "average 
target  system  sighting  efficiency."  This  average  will  converge  to  the  true 
mean  (for  each  point)  as  OPTIMIZE  calls  SIDEN  iteratively. 

Development  of  a  test  procedure  which  could  be  used  by  OPTIMIZE  to 
determine  whether  enough  runs  of  SIDEN  or  BATTLE  have  occurred  is  common 
to  both  cases.  The  common  concern  is  that  the  respective  means  are  within 
specified  intervals  with  a  given  confidence.  Since  it  could  be  expected 
that  the  number  of  iterations  needed  would  be  less  than  30,  the  objective 
is  to  obtain  the  mean  within  +  t j. ,  s/z*  or 
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P(  4"  t*/Zjr\'\  5/y7T<  jX  <  h  +  tj:/2irt-l  3 /-VrT* )  -  I  -  aC 

The  underlying  assumption  is  that  y  is  distributed  normally,  and  this 
is  addressed  below. 

If  the  confidence  interval  within  which  the  true  mean  may  be  found 
90%  of  the  time  is  expressed  in  terms  of  the  sample  mean  as: 

4-  J  4  t  /a  £ 

we  are  saying  that  the  true  mean  is  within  +  10%  of  y  with  a  confidence 
of  .9.  Therefore, 

J  4  =  s/vn 

Vie  also  know  that  (Ref  15284)  an  unbiased  estimator  for  the  sample  variance 

is  given  by  n 

S2  =  — ! —  T  ( Mi  -  hT 
n-»  f- 

t--* 

Therefore,  we  may  say  that 

Increasing  n  from  2  until  the  left  hand  side  is  greater  than  or  equal  to 
the  right  hand  side  provides  the  ,9  confidence  that  the  true  mean  is  in 
the  interval  y  +  .1  y  . 

This  "stopping  rule"  test  is  somewhat  similar  to  that  formulated  by 
Stein  (Ref  17:  479  -  481),  and  can  be  termed  a  "sequential  approach"  to 
determining  the  minimum  number  of  iterations. 

In  the  case  of  the  SUMEFF  values  being  discussed,  the  "target  system 
sighting  efficiency"  values  reflect  averages  over  all  targets,  and  there¬ 
fore  would  tend  to  be  normally  distributed  in  accordance  with  the  central 
limit  theorem  (Ref 17:255),  even  for  averages  over  14  or  fewer  targets. 

If  the  above  defined  "stopping  rule"  were  to  be  used  to  build  some 
confidence  in  the  values  in  the  SUMEFF  matrix,  storage  of  the  data 

needed  to  update  the  sample  mean  and  variance  with  each  call  of  SIDEN 
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would  prove  unmanageable.  Random  sampling  techniques  could  be  used 
to  test  fewer  than  8250  points  with  only  a  slightly  lowered  confidence 
in  the  interval  of  the  mean.  For  example,  random  test  of  about  740  of 
the  8250  elements  of  the  SUMEFF  matrix  would  yield  a  probability  of  .99 
that  at  least  one  of  the  elements  which  has  a  sample  variance  in  the  top 
(highest)  5%  of  all  8250  elements  would  be  found  (Refl6:63).  The  cal¬ 
culation  for  the  required  number  of  random  samples  is  given  by: 

C.z  in  C  I  ”  P  603/  In  (  j~  -f ) 
where  n  =  Number  of  random  samples  (out  of  1000), 

f  *  Fraction  (percent  of  population  of  1000  within  which 
at  least  one  sample  will  be  found  with  a  probability 
of  p  (f). 

In  this  case, 

n*  In  (  I-/M)  /  In  (l-.O.f)  sSl.  7<S 

This  formulation  is  constructed  for  a  population  of  1000,  so  we  want 

n*  C8Z50 / iqoo)  69/75  >  '7*H 

This  procedure  is  the  one  recommended  for  use  with  the  above  defined 
stopping  rule  in  providing  some  confidence  or  reliability  in  the  "average 
target  sighting  efficiency"  values  in  the  matrix  SUMEFF. 

Since  the  developed  stopping  rule  is  based  on  the  assumptions  of 
a  normal  (or"near-normal")  population,  a  "goodness-of-fit"  test  should 
be  run  on  the  output  of  SIDEN  to  provide  some  degree  of  confidence  that 
normality  does  in  fact  exist, 

Pritsker  and  Pegden  ( Re r 1 8:465-470)  discuss  a  "derived  observation" 
which  may  assume  either  a  "0"  or  a  "1"  value.  The  results  of  a  BATTLE 
(variable  "SUCESS"  assumes  a  value  of  1  for  a  "successful"  BATTLE,  and 
a  value  of  0  fo»*  an  "unsuccessful"  BATTLE)  fits  this  concept  very  well. 
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Although  the  stopping  rule  derived  above  could  also  be  applied  here, 
it  would  be  reasonable  to  expect  that  the  sample  variance  of  the 
SUCESS  variable  would  dampen  to  a  small  value  more  quickly  than  some 
of  the  "average  target  system  sighting  efficiencies"  stored  in  SUMEFF. 

If  applicable  tests  should  fail  to  reject  this  hypothesis,  the  stopping 
rule  applied  to  the  SIDEN  outputs  averaged  into  the  SUMEFF  matrix  would 
also  provide  assurance  that  enough  BATTLES  have  been  run  to  provide 
confidence  in  the  success  or  failure  of  the  current  weapon  system  to 
meet  the  mission  requirements  Crto .  of  battles  =  no.  of  SIDEN  calls). 

For  the  purpose  of  expedience,  SIDEN  and  BATTLE  will  each  be  called 
three  times  during  simulation  runs  made  for  model  demonstrations.  The 
results  of  these  experiments  must  therefore  be  evaluated  with  this  in 
mind. 

The  first  area  to  be  investigated  in  this  section  is  the  HELBASE 
output  of  a  weapon  system  which  meets  all  mission  requirements.  The 
output  consists  of  weapons  of  specified  weapon  types  placed  on  specific 
latitudes  or  in  specific  orbits.  The  primary  output  parameter  of  interest 
is  the  total  number  (of  all  weapon  types)  of  weapons  needed  to  complete 
the  mission. 

The  confidence  intervals  imposed  on  the  means  of  SUMEFF  values  and 
the  BATTLE  variable  SUCESS  should  tend  to  make  the  output  weapon  system 
deterministic.  However,  in  cases  where  the  placement  of  the  first 
weapon  involves  the  choice  by  OPTIMIZE  from  several  weapon  paths  with 
nearly  equal  path  efficiencies,  the  overall  number  of  weapons  deployed 
may  change  depending  on  the  placement  of  the  first  weapon. 

This  completes  both  the  stochastic  event  testing  discussion  and 
the  verification  section.  The  two  concepts  of  model  validity  left  are 
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model  credibility  and  confidence  development. 

Credibility 

Rather  than  eliciting  responses  from  experts  by  scenario  description 
and  presentation  of  the  simulation  results,  face  validation  of  HELBASE 
was  accomplished  by  elicitation  of  responses  from  experts  concerning 
the  physical  processes  modeled,  the  attributes  selected,  and  the  design 
structure  of  the  model. 

Several  people,  each  expert  in  a  field  associated  with  the  physical 
processes,  attributes,  or  structure  of  HELBASE,  were  consuUed  in  the 
areas  listed  in  table  IV  .  The  refinements  and  changes  listed  were  made 
at  the  suggestion  of  the  corresponding  expert. 

A  second  aspect  of  the  credibility  part  of  the  validation  process 
is  the  issue  of  sensitivity  analysis.  Model  demonstration  (Chapter  V) 
incorporates  an  experimental  design  which  investigates  the  sensitivity 
of  HELBASE  output  to  changes  in  target  system  negation  time  and  target 
type  negation  percentages. 

An  apropos  technique  for  investigating  the  significance  of  each  of  the 

1/ 

model  inputs  would  be  a  2  factorial  experimental  design  based  on  factors 
which  are  groups  of  the  22  individual  variables.  Even  though  this  would 
be  a  time-consuming  process,  it  would  definitely  have  value  for  the 
serious  model  user. 

Confidence 

A  useful  tool  in  building  confidence  in  the  performance  of  submodels 
PROP  and  EFFECT  is  comparison  of  "one-on-one"  (one  weapon  firing  at  one 
target)  results  with  another  model.  Such  models  exist  in  varying  degrees 
of  complexity,  one  fairly  accurate  one  built  by  Peckham  and  Davis  of  the 
Air  Force  Weapons  Laboratory  (AFWL)  (Ref 20).  For  the  "first  order 
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TABLE  IV 

Face  Validation  References 

Dunne,  Edward  J. ,  Lt  Col  —  Model  structure.  Sighting  efficiency  attribute. 

Havcy  »  Janies  H. ,  Lt  Col  —  Physical  process  of  propagation.  Definition  of 

weapon  parameters. 

Torvik,  Peter  J.  —  Laser  effects  physical  process. 

\ 

Wiesel,  William  E. ,  Capt  —  Orbital  relationships. 
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accuracy"  of  HELBASE,  an  accuracy  of  plus  or  minus  25%  of  the  AFWL 
model  output  would  seem  reasonable.  However,  this  determination  must 
be  made  by  the  user  depending  on  the  nature  of  his  decision  and 
experimental  design.  If  necessary,  adjustments  to  submodels  PROP  and 
EFFECT  could  be  made. 

Finally,  the  potential  user  will  find  that  documentation  of 
HELBASE  is  fairly  thourough.  A  "documentation  list"  provided  by 
Shannon  (Ref  19:  262)  calls  for  the  following  items. 

Flow  diagrams  of  each  module  and  the  overall  model  — 
Chapter  II. 

Description  of  inputs  necessary  for  executing  the  program, 
including:  input  card  number,  symbol,  definition, 
whether  integer  or  real,  and  the  field— Chapter  V. 

Definition  of  program  variables  not  used  as  inputs  to  the 
program--  Appendix  C. 

Verbal  descriptions  of  all  modules  as  to  purpose  and 
function—  Chapter  II. 

Input  deck  setup  to  run  on  the  computer  of  interest- 
provided  separately. 

Listing  of  the  program--  Appendix  C. 


t 
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V.  Use  of  the  Model 

The  purpose  of  this  chapter  is  to  demonstrate  the  use  of  HELBASE 
tn  two  fairly  typical  applications,  and  to  discuss  the  use  of  the  model 
t  in  any  proper  experiment.  The  general  use  of  the  model  is  covered  in 

the  format  of  an  informal  "user's  guide". 

«  HELBASE  Demonstration 

Two  experimental  designs  were  devised  for  the  model.  The  first 
experiment  demonstrates  the  use  of  HELBASE  to  pick  an  "optimum"  orbital 
radius  for  a  specific  weapon.  The  second  experiment  is  a  sensitivity  analysis 
concerning  two  of  the  HELBASE  input  variables,  an  so  explores  user  concerns 
for  a  specific  model  application. 

Experiment  #1  --  One-Dimensional  Optimization  of  an_  Orbital  Radius. 

Of  the  22  variables  which  may  be  input  by  tne  user,  21  are  held  constant 
for  this  experiment.  The  only  variable  is  the  orbital  altitude  of  the 
space-based  weapon  type  being  employed. 

The  hypothetical  decision  here  is  that  of  choosing  the  best  orbital 
radius  for  the  < pace-based  weapon  type.  An  in-place  ground-based  weapon 
type  is  deployed  as  shown  in  Table  V  These  ground-based  weapons 
constitute  an  in-being  force,  and  will  be  used  in  conjunction  with  the 
space-based  weapon  type  in  performing  the  mission.  The  continued  exist¬ 
ence  and  use  of  these  ground-based  weapons  is  certain. 

*  Target  type  #1,  which  is  6  targets  in  250  NM  slightly  elliptical 

orbits,  was  chosen  as  the  target  system.'  A  target  type  negation  percent¬ 
age  of  .,66  translatesto  4  out  of  the  6  targets  being  negated  for  mission 
denial.  A  6?)  minute  target  system  negation  time  gives  the  in-place 
ground-based  system  of  2  weapons  and  the  space-based  weapons  just  1  hour 
to  negate  the  4  targets.  Sorting  of  targets  of  opportunity  by  least 
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Table  V 

User  Inputs  for  Experiment  #1 


Ground  Based  Weapon  Type  #1 


Wavelength  (OF) 

HEL  Power  at  Source 

HEL  Waist  Size  at  Source 

Maximum  Firing  Cycles  Allowed 

Maximum  Firing  Time  per  Operating  Cycle 

Recycle  Time 

Maximum  Range  Allowed 

Probability  of  Failure 

Mode  Number 

Altitude 

Absorption  Coefficient 


3800  NM 
250,000  Watts 
3  Meters 
20 

10  Seconds 
120  Seconds 
c*>  Default 

0 

1 

0  -1 
.07  KM  1 


Wavelength 

HEL  Power  at  Source 

HEL  Waist  Size  at  Source 

Maximum  Firing  Cycles  Allowed 

Maximum  Firing  Time  per  Operating  Cycle 

Recycle  Time 

Maximum  Range  Allowed 

Probability  of  Failure 

Mode  Number 

A1 ti tude 

Absorption  Coefficient 


10,600  NM 
50,000  Watts 
1  Meter 
5 

20  Seconds 
60  Seconds 
o» 


500  NM  -  (17,812  KM) 

.143  KM"1 


#1  (6  targets  in  250  NM  slightly  elliptical  orbits) 


Mission  Data 
Target  Type 

Target  Type  Negation  Percentage 
Target  System  Negation  Time 
Target  Type  Priority 


Battle  Management  Data 


Target  Type 

Firing  Priorities 

Sorting  by  Irradiation  Time 


1 

.66  (4  targets) 
60  minutes 
Default  (1) 


1 

Default  (4) 
Yes  (1) 


irradiation  time  is  employed  as  a  tactic.  Obviously,  assigning  firing 
priorities  is  not  appropriate  since  only  one  target  type  was  selected 
for  this  experiment. 

The  decision  makers  are  faced  only  with  the  selection  of  the  circular 
orbit  radius  of  the  space-based  weapons  in  this  case.  The  decision  has 
already  been  made  to  deploy  the  space  weapons  to  complement  the  in-being 
ground  weapons,  and  all  of  the  space  weapons  are  to  be  placed  at  the  same 
orbital  radius.  The  objective  is  to  find  the  orbital  radius  for  the  space 
weapons  which  causes  the  fewest  number  of  space  weapons  to  be  needed  to 
meet  the  mission  requirements. 

The  experimental  design  is  based  on  the  one-dimensional  Golden-Section 
search  technique  (Ref  16:32-35).  Of  the  class  of  one-dimensional  search 
techniques  based  on  the  assumption  of  unimodality,  the  Golden  Section 

technique  is  the  second  most  efficient  in  terms  of  the  number  of  obser- 

.  \ 

vations  needed.  It  is  only  slightly  less  efficient  than  the  Fibonacci 
technique.  The  observations  consist  of  the  number  of  space  weapons  re¬ 
quired  along  with  the  radius  of  the  weapons  input.  So  the  "x"  (independ¬ 
ent)  variable  is  the  input  radius,  and  the  "f(x)"  (output)  response  is 
the  number  of  space  weapons  required  to  meet  the  mission  requirements. 

As  stated  above,  the  Golden  Section  search  technique  is  based  on  the 
assumption  that  the  response  "f(x)"  is  unimodal  over  the  range  of  "x" 
being  considered.  In  other  words,  the  f(x)  response  must  exhibit  one 
of  four  kinds  of  behavior  as  x  varies  between  its  minimum  and  maximum 
values.  The  response  f(x)  must:  1)  strictly  increase,  2)  strictly 
decrease,  3)  strictly  increase  to  a  maximum  and  then  strictly  decrease, 
or  4)  strictly  decrease  to  a  minimum  and  then  strictly  increase.  In  all 
but  the  third  case,  the  optimum  solution  will  be  at  the  minimum  or  maxi¬ 
mum  value  of  x. 
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The  literature  does  not  offer  a  method  of  estimating  the  unimodality 
of  a  simulation  response  of  an  unknown  character.  There  is  a  rough 
intuitive  reasonableness  in  assuming  unimodality  in  this  case.  It  seems 
justifiable  that  one  and  only  one  maximum  would  exist  for  the  HELBASE 
response  for  a  range  of  X  from  500  to  17,872  KM.  If,  on  the  other  hand, 
the  target  type  chosen  had  been  type  #3  (500  NM,  926.5  KM  circular  orbits), 
and  X  were  allowed  to  vary  from  500  to  17,872  KM,  an  assumption  of  uni¬ 
modality  would  not  seem  reasonable.  The  existence  of  local  minima  or 
maxima  as  X  moves  from  below  to  above  the  target  system  altitude  would 
seem  probable. 

In  the  case  under  consideration,  the  target  system  is  at  463.25  KM, 
and  the  weapon  altitude  is  allowed  to  vary  from  500  to  17872  KM.  So 
local  minima  or  maxima  would  not  be  expected.  As  a  rudimentary  check 
on  the  unimodality  assumption,  the  response  of  HELBASE  can  be  inspected 
during  the  Golden  Section  search  procedure. 

This  procedure  is  based  on  the  simple  premise  of  iteratively  de¬ 
creasing  the  size  of  the  feasible  interial  of  X  until  an  optimum  value 
X*  can  be  approximated  with  some  degree  of  accuracy.  The  initial 
inter  al  (500,  17872)  is  divided  into  overlappina  subintervals  of  length 
.618  of  the  interval  length.  The  two  subintervals  are  measured  from 
opposite  ends  of  the  initial  interval.  That  is,  point  1  is  located 

10735.9  KM  from  the  lower  end  of  the  interval,  and  point  2  is  located 

10735.9  KM  from  the  upper  end  of  the  interval.  For  the  initial  interval, 

point  1  is  at  11235.9  KM,  and  point  2  at  7136.1  KM.  HELBASE  responses 

(number  of  weapons)  are  then  obtained  for  each  altitude.  If  the  HELBASE 

response  for  an  orbital  radius  of  71 36 . 1  KM  ( f ( 71 36 . 1 )  were  less  than 

f ( 11235.9) ,  the  region  of  the  interval  (500,  17872)  from  500  to  7136.1 
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need  not  be  considered  further.  So,  we  take  as  the  next,  inter  al  to 
be  evaluated  (7136.1,  17872).  Similarly,  if  f( 71 36.1)  were  greater  than 
f ( 1 1235. 9) ,  the  region  of  (500,  17872)  from  11235.9  to  17872  need  not 
be  considered  further,  and  we  would  take  as  the  next  interval  (500), 
11235.9).  Note  that  in  both  cases,  the  initial  interval  would  be 
decreased  38.2%  (.618  of  the  initial  length),  and  that  one  of  the  res¬ 
ponses  needed  for  the  evaluation  of  the  next  interval  would  have  already 
been  obtained. 

The  number  of  iterations  of  the  above  procedure  is  related  to  the 
size  of  the  "final"  interval  (within  which  the  final  value  of  X  will 
He)  is 

N -  /”  ^  ■  j 

in  *  C,iS  ,  where 

N  is  the  number  of  model  responses  needed, 

Ln  is  the  length  of  the  final  (Nth) 
interval ,  and 

Lo  is  the  length  of  the  initial  interval. 

If  the  object-’ ve  is  to  find  a  final  interval  of  length  1000  KM  which  is 
guaranteed  to  contain  the  maximum  X*  under  the  assumption  of  unimodality, 

6  model  responses  will  be  needed.  That  is,  /,»  (/oco> -in  (ni72) 

- - - - s' 

in  (•  6  iS) 

In  addition,  a  refinement  of  the  final  X  aoproximation  to  X*  can  be  found 
by  graphing  the  functional  responses  against  the  input  X's  and  interpolating 
(or  extrapolating)  through  the  final  interval. 
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Analysis  of  the  results  of  this  experiment  is  difficult  due  to 
the  single  variable  nature  and  the  small  number  of  responses.  The 
"inference  space"  of  concern  here  is  fairly  restricted.  The  "optional 
radius"  found  suffers  several  degrees  of  approximation.  It  endures 
not  only  the  degree  of  approximation  caused  by  HELBASE  assumption, 
accuracy,  scope  and  limitations,  but  an  additional  degree  inferred  by 
the  unimodality  assumption.  Part  of  the  implied  inference  space  for 
the  results  of  this  experiment  is  that  the  optional  radius  found  can 
only  be  "trusted"  for  the  system  configuration  used.  That  is,  the 
number  and  placement  of  the  space  weapons  cannot  be  assumed  to  be 
independent  of  the  number,  characteristics  and  placement  of  the  ground 
weapons.  Simple  tests,  such  as  statistical  significance  tests  of  vary¬ 
ing  the  ground-based  weapon  parameters  could  be  used  to  reject  the 
significance  of  changes  in  those  parameters.  In  this  way  an  extension 
of  the  inference  space  in  which  the  optimum  radius  result  could  be 
considered  valid  could  be  obtained. 

Experiment  #2  —  Two  Factor  Sensitivity  Analysis.  In  general ,  the 
interrelationship  between  any  two  factors  in  HELBASE  is  dependent  upon 
the  values  assigned  to  the  other  20  factors.  Therefore,  sensitivity 
analysis  can  only  have  meaning  within  defined  ranges  of  the  factors  of 
interest  with  the  rest  of  the  input  variables  being  held  constant. 

The  two  factors  being  considered  here  are  the  target  system  negation 
time  and  the  target  type  negation  percentage.  In  the  decision  situation 
described  below,  the  target  system  negation  time  will  be  varied  from 
1  hour  and  30  minutes  to  two  hours,  and  the  negation  percentage  for  both 
target  types  will  be  varied  from  .60  to  .70, 


The  mission  envisioned  in  this  case  is  to  deny  a  potential  enemy 
the  use  of  a  target  system  consisting  of  3  satellites  in  50(5  NM  (926.5 
KM)  circular  orbits  (type  3)  and  6  satellites  in  250  NM  (463.25  KM) 
slightly  elliptical  orbits  (type  1).  The  only  HEL  weapon  available  for 
deployment  is  a  deuterium  flouride  system  nearing  operational  testing  as 
a  ground  battle  station. 

Ihe  decision  has  been  made  to  deploy  as  many  of  these  systems  as 
are  "needed". 

Weapon  parameters  are  listed  in  Table  VI. 

It  has  beer,  decided  that  it  is  roughly  "10  times  more  imoortant" 
to  negate  the  type  1  targets  than  the  type  2.  Therefore,  a  target  type 
priority  of  10  has  been  assigned  to  the  6  type  1  targets.  In  conjunction 
with  this,  target  type  1  is  assigned  a  firing  priority  of  1.  The  tactic 
of  firing  at  a  target  which  requires  a  shorter  irradiation  time  for  tar¬ 
get  negation  is  used. 

The  only  decision  variables  not  "nailed  down"  at  this  point  are  the 
negation  percentage  and  the  target  system  negation  time. 

The  decision  situation  is  as  follows.  We  will  assume  that  HELBASE 
was  run  with  the  target  type  negation  percentage  set  to  .7  and  the  target 
system  negation  time  set  to  1  hour  and  30  minutes.  The  model  output 
indicated  that  x  weapons  were  needed  on  various  latitudes.  A  review 
of  this  deployment  scheme  yielded  a  concern  over  whether  the  x  number  of 
weapons  could  be  reduced  by  increasing  the  allowed  mission  time  or  by 
decreasing  the  percentage  of  each  target  type  required  to  be  negated,  or 
perhaps  some  of  both. 

The  lowest  percentage  that  the  decision  maker  will  accept  as  reason¬ 
able  is  .6,  and  the  maximum  target  system  negation  time  is  2  hrs.  A  .6 
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Table  VI 


User  Inputs  for  Experiment  #2 


DF  Weapon 


Wavelength 

HEL  Power  at  Source 

HEL  Waist  Size  at  Source 

Maximum  Firing  Cycles  Allowed 

Maximum  Firing  Time  per  Operating 

Recycle  Time 

Maximum  Range  Allowed 

Probability  of  Failure 

Mode  Number 

Altitude 

Absorption  Coefficient 


3800  NM 

50,000,000  watts 
3  meters 
10 

Cycle  10  seconds 

60  seconds 
5000  KM 
0 
1 

.0? 


Target  Types  #1  and  #3. 
Mi ssion  Data 


Target  Type  #1 : 

Target  Type  Negation  Percentage 

.60  to  .70 

Target  System  Negation  Time 

90  to  120 

Target  Type  Priority 

10 

Target  Type  #3: 

Target  Type  Negation  Percentage 

.60  tc  .70 

Target  System  Negation  Time 

90  to  120  min. 

Target  Type  Priority 

1 

Battle  Management  Data 

Target  Type  #1 : 

Firing  Priority 

1 

Target  Type  #2 : 

Firing  Priority 

4 

Sorting  by  irradiation  Time:  Yes 
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percentage  translates  to  2  out  of  the  3  type  3  satellites, 

(  .6  X  3  -  l.a  a  t  \  and 

4  out  of  the  6  type  1  satellites.  The  .7  percentage  translates  to  all  of 
the  type  3  satellites  and  S  of  the  6  type  1  satellites.  An  examination 
of  these  alternatives  reveals  that  we  are  actually  concerned  with  .66 
percentages  for  the  low  end  and  .83  (type  1)  and  1.00  (tyoe  3)  for  the 
high  end.  This  small  revelation  is  even  more  palatable  to  the  decision 
maker,  since  he  does  not  have  to  accept  an  actual  negation  percentage  of 
.6,  and  the  target  type  negation  percentage  alternatives  are  set  at  .6 
and  .7  for  both  target  types. 

The  decision  maker  regards  the  increasing  of  mission  time  to  two  hours 
and  the  decreasing  of  negation  percentage  to  .6  as  equally  degrading  to 
the  mission.  Therefore,  an  "equi-valued"  pair  of  alternatives  is  as 
follows: 

Alternative  #1:  Negation  percentage  at  .6, 
mission  time  at  90  minutes. 

Alternative  #2:  Negation  percentage  at  .7 

mission  time  at  120  minutes. 

Since  we  have  two  factors  at  two  levels,  it  would  seem  reasonable  to 
pursue  a  factorial  experimental  design.  A  full  factorial  two  factor- 
two  level  design  (Ref  B  =  163-165)  requires  that  we  also  evaluate  the  "high- 
high"  and  "lew-low"  combination  in  order  to  allow  assessment  of  factor  signi 
ficance  and  any  interrelationshios,  between  the  factors.  Therefore,  we  need 
two  more  experimental  alternatives: 

Alternative  #3:  Negation  percentage  at  .7, 
mission  time  at  90  minutes, 

Alternative  #4:  Negation  percentage  at  .6, 
mission  time  at  120  minutes. 
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Of  course,  alternative  #3  is  the  original  factor  combination  which 
produced  the  X  weapon  requirement.  This  data  point  can  be  used  as  one 
of  those  needed  for  this  experiment.  Inclusion  of  these  alternatives  in 
the  full  factorial  design  will  allow  an  evaluation  of  the  other  stated 
alternative  "or  perhaps  some  of  both". 

The  experimental  design  is  depicted  in  figure  21  .  Since  the  HELBASE 
output  (number  of  weapons)  would  be  expected  to  yield  fairly  constant 
responses  for  any  combination  of  the  factor  levels,  Three  responses 
for  each  combination  will  be  obtained.  Therefore,  12  model  runs  will 
be  needed  in  order  to  evaluate  all  combinations  of  factors  with  two  repli¬ 
cations.  See  Table  VII  for  a  listing  of  the  input  cards. 

2 

The  implied  model  underlying  this  2  full  factorial  experimental  design 
is 

-  JJ.-*  o£‘t  +  fi.  Gijk 

where 

i  is  the  level  of  the  negation  percentage, 
j  is  the  level  of  mission  time, 
k  is  the  kth  observation  at  factor  levels  T,  3, 
y*is  the  mean  of  all  responses, 

<<'is  the  contribution  of  negation  percentage 
at  level  i, 

is  the  contribution  of  mission  time  at  level  j. 

Is  the  contribution  of  the  interaction  between 
negation  percentage  at  level  i  and  mission  time 
at  level  j,  and 

is  a  random  error  term  (unexplained  contribution). 

The  analysis  of  the  experimental  results  is  done  in  a  analysis  of 
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TARGET  SYSTEM  NEGATION  TIME 


90  ito 


TARGET  .6 

TYPE 

NEGATION 

xin 

X1 1 2 

X113 

x121 

*122 

X123 

PERCENTAGE 

x211 

X221 

.7 

*212 

X222 

x213 

x223 

Source  of  Variation  Degrees  of 
"  "  ”  Freedom 


Negation  Percentage  1 

Mission  Time  1 

Interaction  1 

Error  8. 

Total  11 


Sum  of 
Squares 

Mean 

Squares 

F 

SSA 

msa=ssa 

MSA/MSE 

SSB 

msb=ssb 

MSB/MSE 

ssAb 

msab=ssab 

MSAB/MSE 

SSg 

MSE=SSE 

a 

Fig  21.  Experiment  #2  Design 


TABLE  VII 

Input  Cards  Used  for  Experiment  #2 

Note:  All  values  are  real,  and  all 
cards  are  free-formatted. 


LOW-LOW 


CARP  #  FIRST  REPLICATION  (LOW- LOW) 

1  8972. 

2  1. 

3  3800  50000000.  3.  10.  10.  60.  5000.  0.  1.  0.  .07 

4  1. 

5  3. 

6  99. 

7  1.  .6  90.  10. 

8  3.  .«  90.  1. 

9  1.  1.  1. 

10  3.  4.  1. 

11  7/8/9 

Legend  for  the  Above  Data  Set 

CARD  %  LEGEND 

1  (8972.)  Random  number  seed. 

2  (1.)  Number  of  weapon  tyoes  input. 

3  (3800.)  Wavelength  ’/i  nanometers. 

(50000000.)  Weapon  power  in  watts. 

(3.)  HEL  waist  size  in  meters. 

(10.)  Maximum  allowed  firing  cycles. 

(10.)  Maximum  firing  time  per  cycle  in  seconds. 

(60.)  Recycle  time  in  seconds. 

(5000.)  Maximum  range  in  kilometers. 

(0.)  Probability  of  weapon  failure. 

(1.)  Mode  number. 

(0.)  Weapon  altitude  in  kilometers. 

(.07)  Absorption  coefficient  in  kilometers"  . 

4  (1.)  Target  type  1  selected. 

5  (3.)  Target  type  3  selected. 

6  (99.)  Target  type  selection  terminator. 

7  (1.)  Target  type  1. 

(.6)  Negation  percentage  for  target  type  1. 

(90.)  Target  system  negation  time  in  minutes. 

(10.)  Target  ype  1  priority. 

8  Same  as  car^  #7. 

9  (1 . )  Target  type  1 . 

(1.)  Firing  priority  against  target  type  1. 

( 1 . )  Sorting  by  irradiation  time  selected. 

10  Same  as  card  #9. 

11  End  of  record. 
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CARD  # 

SECOND  REPLICATION 

12 

5555. 

13 

• 

• 

Same  as  above. 

22 

CARD  # 

third  Replication 

23 

14670927. 

24 

• 

• 

Same  as  above. 

33 

LOW- HIGH 

CARD  # 

first  replication 

34 

4492. 

35 

1. 

36 

3800.  500001 

37 

1. 

38 

3. 

39 

99. 

40 

1.  .6  120.  ' 

41 

3.  .6  120. 

42 

1.  1.  1. 

43 

3.  4.  1. 

44 

7/8/9 

CARD  # 

SECOND  REPLICATION 

45 

92929292. 

46 

• 

• 

Same  as  above. 

55 

CARD  # 

THIRD  REPLICATION 

56 

197555. 

57 

• 

Same  as  above. 

(LOW-LOW) 


:  LOW-LOW) 


(LOW-HIGH) 


5000000(3.  3.  10.  10.  60.  5000.  0.  1.  0..0  7 


(LOW-HIGH) 


:  LOW-HIGH) 
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SECOND  REPLI CAT I ON 
7777. 

Same  as  above. 


THIRD  REPLICATION 
9876543. 


Same  as  above. 


variance  format,  discussed  widely  in  the  literature  (Refs  15:  472- 
477;  19:  162). 

User’s  Guide 

The  above  two  experimental  designs  serve  to  exphasize  some  aspects 
of  using  the  HELBASE  model  in  the  support  of  specific  decisions. 

Because  HELBASE  was  conceived  and  designed  to  be  a  flexible  tool 

for  the  decision  maker,  analyst  or  manager,  amenable  to  a  variety  of 

>‘  • 

experimental  designs  in  support  of  a  wide  range  of  decisions,  a  certain 
amount  of  responsibility  for  proper  use  lies  with  the  user. 

In  addition  to  the  normal  precautions  and  care  which  should  be  taken 
with  the  use  of  any  model,  the  HELBASE  user  must  be  aware  of  certain 
pitfalls  and  biases  which  must  be  warded  against.  These  biases  can 
occur  due  to  improper  use  of  the  22  input  variables,  or  through  mis¬ 
understanding  of  the  inference  space  surrounding  the  model  output. 
Accordingly,  this  informal  "user's  guide"  is  divided  into  two  sections-- 
"Pre-Analysis  and  Experimental  Design",  and  "Post-Analysis". 

Pre-Analysis  and  Experimental  Design.  The  experimental  designs 
available  to  the  user  are  not  limited  by  the  HELBASE  model  structure. 
However,  after  the  user  thoroughly  understands  the  decision  situation 
far  which  the  experiment  will  be  carried  out,  has  satisfactorily  expressed 
the  experimental  objective,  and  has  designed  the  experiment,  a  few  other 
considerations  should  be  made  before  the  experiment  is  run.  These 
considerations  are  related  to  the  physical  processes  modeled,  the  inter¬ 
actions  among  the  input  variables,  and  the  sensitivity  of  the  model 
output  to  the  decision  variables  of  concern. 

First,  the  user  must  have  confidence  in  PROP.  Especially  critical  is 
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a  check  on  the  possibility  of  thermal  blooming  for  ground-  or  air-based 
weapons.  In  addition  to  the  thermal  blooming  approximation  technique 
referenced  in  Chapter  I,  a  propagation  model  such  as  the  one  proffered 
by  Peckham  and  Davis  ( Ref  2C$  can  be  used  to  build  confidence  in  the 
HELBASE  submodel  PROP  in  any  specific  application. 

An  additional  check  in  this  area  would  be  assurance  by  the  user  that 
the  lasc  ‘  power,  wavelength  and  location  combinations  input  will  yield 
a  lethal  effect  on  a  significant  number  of  the  targets.  HELBASE  cannot 
choose  a  "most  efficient"  path  if  no  weapon  negations  are  feasible.  In 
a  case  such  as  this,  the  model  would  terminate  execution  when  OPTIMIZE 
attempted  to  find  the  maximum  of  a  set  weapon  path  efficiencies  which  are 
all  equal  to  zero. 

Second,  the  values  of  all  non-decision  variables  in  the  22  input 
variable  set  must  be  chosen  with  care.  Many,  if  not  most,  can  be  assign¬ 
ed  values  which  are  reasonable  for  the  decision  situation  under  consider¬ 
ation. 

Due  to  the  generally  unexamined  nature  of  interrelationships  among 

v 

the  22  input  variables  with  respect  to  the  HELBASE  output  weapon  system, 
even  acceptance  of  model  default  values  should  be  done  with  care.  For 
example,  acceptance  of  the  default  target  type  negation  percentage  of  .8 
could  very  well  affect  the  significant  sensitivity  range  of  the  target 
system  negation  time.  If  the  sensitivity  of  the  HELBASE  output  to  the 
mission  time  were  inportant  to  the  user's  experiment,  this  unwanted  and 
unknown  (to  the  user)  bias  could  well  affect  model  output  as  well  as  the 
resulting  decision. 

The  general  area  of  sensitivity  analysis  is  the  third  and  last  area 

in  which  some  caution  should  be  shown  by  the  user.  Since  a  generalized 
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sensitivity  analysis  and  definition  of  interrelationships  among  all 
22  input  variables  would  prove  an  intolerable  computational  feat,  the 
potential  user  may  find  it  advisable  to  conduct  a  restricted  sensi¬ 
tivity  analysis  over  the  decision  variables  chosen  for  the  specific 
experiment.  An  approximation  of  the  individual  and  joint  effects  can 
also  De  calculated  as  described  in  experiment  #2. 

PoSt-Analysis.  In  conjunction  with  the  traditional  statistical 
analysis  performed  for  the  user's  experimental  design,  two  additional 
areas  of  concern  should  be  considered  after  the  experiment  has  been  run. 

Because  the  normal  purpose  of  a  simulation  experiment  is  to  aid  a 
decision  maker,  the  user  should  insure  that  the  inference  space  (the 
combination  of  model  assumptions,  limitations,  scope  and  accuracy  with 
the  user's  specific  choice  of  constant  as  well  as  decision  variables 
within  the  experimental  design)  has  been  extended  and  is  valid  over 
the  decision  situation. 

The  second  major  consideration  that  may  have  value  to  a  user  is 
an  analysis  to  selectively  improve  upon  the  "baseline"  HELBASE  weapon 
deployment  scheme.  The  output  of  the  model  is  termed  "baseline"  because 
some  judicious  adjustments  of  the  weapon  locational  parameters  produced 
by  HELBASE  can,  in  many  cases,  improve  the  performance  of  the  weapon 
system.  As  discussed  in  Chapters  II  and  III,  HELBASE  places  ground- 
and  air-based  weapons  on  latitudes,  and  defines  only  five  of  the  six 
orbital  parameters  for  space-based  weapons. 

This  means  that,  for  latitudes  with  mere  than  one  weapon  assigned, 

the  exact  locations  of  the  weapons  around  the  latitudinal  circle  is 

treated  randomly  by  the  model,  with  no  effect  on  the  "baseline"  capability 

of  the  weapon  system  to  satisfy  the  mission  requirements.  For  example, 
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consider  a  HELBASE  output  which  places  three  ground-based  weapons  at 
theta  *  30°  (60  N  latitude).  Random  treatment  of  the  locations  of 
the  three  weapons  on  the  60°  North  latitude  implies  that  the  three 
weapons  may  be  placed  virtually  anywhere  around  the  latitudinal  circle 
while  retaining  the  "baseline"  capability  to  meet  the  mission  require¬ 
ments.  Therefore,  thoughtful  placement  of  the  three  weapons  may  yield 
significant  improvement  in  weapon  system  performance.  In  this  example, 
distributing  the  three  ground-based  weapons  at  120°  intervals  around 
the  latitudinal  circle  could  be  expected  to  improve  the  performance  of 
the  weapon  system. 

Similar  improvements  may  be  expected  for  air-based  weapons.  In 
addition  to  the  distribution  of  weapons  around  the  latitudinal  circle 
described  above,  air-based  weapons  may  be  flown  off  the  latitudinal 
path  so  as  to  more  efficiently  intercept  targets  in  response  to  the 
known  target  locations  at  the  initiation  of  hostilities.  This  technique 
would  be  expected  to  yield  improvement  in  weapon  system  performance. 

The  basic  reason  behind  this  type  of  performance  improvement  is 
that  the  model  output  weapon  system  is  capable  of  satisfying  mission 
requirements  with  the  time  of  hostilities  initiation  being  random. 

The  unspecified  parameter  in  the  HELBASE  selection  of  space-based 
weapons  is  the  true  anomaly  at  epoch.  This  sixth  orbital  parameter  is 
again  treated  randomly  by  HELBASE  during  the  building  of  the  output 
weapon  system.  A  user  may  find  it  advantageous  to  define  true  anomalies 
for  space  weapons  which  cause  "phase  relationships"  among  the  weapons. 

The  effect  on  weapon  system  performance  of  making  the  above  adjust¬ 
ments  in  the  deployment  parameters  may  be  measured  by  inputting  this 
weapon  system  into  a  BATTLE  and  noting  any  improvement  (decrease)  in 
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the  time  taken  to  satisfy  the  mission  requirements.  It  must  be  noted 
that  some  modification  of  the  current  model  would  be  necessary  to  allow 
this  type  of  HELBASE  employment. 

This  chapter  has  examined  two  experimental  designs  which  would  be 
considered  typical  for  use  with  HELBASE.  In  addition,  considerations 
which  could  improve  the  efficiency  and  validity  of  any  use  of  HELBASE 
were  offered  in  the  "User's  Guide"  section.  The  last  chapter  summarizes 
the  thesis  effort,  investigates  the  implications  and  possible  value  of 
the  deployment  of  the  model,  and  offers  recommendations  to  improve  its 

usefulness. 
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VI  Conclusions  and  Recommendations 


The  most  important  consideration  here  is  whether  the  thesis  effort 
met  its  objective,  and  whether  the  HELBASE  model  serves  its  stated  pur¬ 
poses.  In  addition  t  these  points,  the  potential  utility  of  the  model 
is  evaluated. 

Recommendations  for  model  improvement  are  many.  The  most  approp¬ 
riate  are  reviewed  here,  and  a  representative  portion  of  the  rest  is  left 
to  Appendix  B. 

Summary 

From  the  thesis  objective  and  problem  statement,  the  model  HELBASE 
was  conceptualized  in  modular  form.  The  physical  processes  modeled  by 
PROP  and  EFFECT  formed  the  accuracy  lower  bound,  and  the  rest  of  the  model 
was  constructed  with  this  in  mind. 

In  addition  to  these  physical  processes  concerning  laser  beam  prop¬ 
agation  and  the  effect  on  the  target,  three  basic  things  were  needed  to 
meet  the  thesis  objective.  An  appropriate  quantitative  attribute  to  meas¬ 
ure  the  "efficiency"  of  weapon  placement  was  developed  and  is  employed  in 
submodel  SIGHT.  A  methodology  for  using  this  attribute  in  picking  a  "most 
efficient"  weapon  path  was  then  constructed.  This  methodology  is  the  basis 
of  submodels  OFTIMIZE,  SIPEN,  TARGET,  and  SPHERE.  After  the  "most  efficient" 
path  was  selected  and  a  weapon  placed  in  this  path,  a  way  to  measure  the 
degree  to  which  this  weapon  and  all  others  previously  placed  could  satisfy 
mission  requirements  was  needed.  Submodel  BATTLE  was  therefore  developed. 

The  aggregation  of  these  submodels,  along  with  an  input  section,  forms 
the  model  HELBASE.  An  external  program  called  ORB  was  written  for  the  sole 
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purpose  of  defining  orbital  points  (see  Appendix  C), 

Validation  and  verification  of  the  model  was  carried  out  with 
the  exception  of  stochastic  tests  on  the  variables  generated  by  sub¬ 
model  SIDEN  and  the  distribution  of  the  BATTLE  response. 

The  model  was  employed  in  support  of  specific  hypothetical  prob¬ 
lems.  In  this  way,  the  complexities  of  HELBASE  use  were  explored  and 
discussed. 

Thesis  Objective  and  Model  Purpose 

As  stated  in  Chapter  I,  "The  purpose  of  this  thesis  is  to  properly 
define  the  appropriite  attributes  of  merit  and  to  develop  a  simulation 
model  which  would  give  a  decision  maker  aflexible  tool  in  investigating 
the  effects  of  these  attributes  on  specific  space  defensive  missions." 

The  model  HELBASE  is  based  upon  the  efficiency  attribute  developed 
in  Chapter  III.  The  range  of  model  applications  can  be  appreciated  by 
inspection  of  the  uses  shown  and  discussed  in  Chapter  V. 

Model  Utility 

The  potential  utility  of  the  model  must  be  judged  in  comparison 
with  other  existing  tools  which  could  be  used  for  the  same  purpose.  A 
comprehensive  literature  search  yielded  no  "many-on-many"  models  of  this 
nature.  It  would  therefore  have  to  be  concluded  that  HELBASE  can  only 
add  to  the  overall  understanding  of  the  "system"  being  investigated. 
Better  understanding  of  the  long  range  impact  of  design  and  development 
decisions,  deployment  options  and  employment  tactics  seem  clearly  justi¬ 
fiable.  In  fact,  with  proper  pre-  and  post-analysis  (See  Chapter  V),  the 
number  of  potential  uses  of  HELBASE  seems  almost  unlimited. 


Recommc  relations 


The  recommenutions  discussed  here  will  be  limited.  See  Appendix  B 
for  an  expanded  discussion. 

Target  Types.  Obviously  the  use  of  HELBASE  by  a  serious  analyst  would 
require  the  changing  of  the  target  types  from  the  four  representative  cases 
to  "real"  target  types  selected  by  the  user. 

Space-Based  Weapons  Orbital  Selection.  Allowing  the  input  of  specific 
elliptical  orbits  for  weapons  would  seem  to  have  value.  Expansion  of  the 
generalized  optimization  technique  to  include  elliptical  orbits  is  discussed 
in  Appendix  B. 

Battle  Tactics.  Additions  to  the  available  options  in  battle  manage¬ 
ment  data  would  give  a  user  more  tactical  options.  Appendix  B  includes 
discussions  of  an  option  to  reserve  battle  stations  to  be  committed  at  a 
specific  time  after  the  initiation  of  hostilities,  and  an  option  to  prioritize 
the  commitment  of  available  battle  stations. 

Cost  Considerations.  In  addition  to  a  look  at  accounting  for  the  fixed 
costs  associated  with  deploying  the  first  weapon  of  a  given  weapon  type,  an 
alternate  "sighting  efficiency"  attribute  is  discussed.  See  Appendix  B. 


108 


Bibliography 


1.  "Technical  Survey:  Particle  Beams,  Laser  Weapons  —  IV  Aviation 
Week  and  Space  Technology,  July  28,  1980. 

2.  Dunne,  Edward  J.  and  Dutton,  Jeffrey  L.  "A  Tool  for  the  Study  of  High 
Energy  Lasers  in  a  Space  Defense  Role".  School  of  Engineering,  Air 
Force  Institute  of  Technology,  Wright-Patterson  Air  Force  Base, 

October  1980. 

3.  Frye,  Alton.  "U.  S.  Soace  Policy:  An  Example  of  Political  Analysis", 
from  Systems  Analysis  and  Pol  icy  Planning:  Applications  in  Defense. 

New  York:  El sevier  North-H'o  lland,  Inc.,  1 9 77. 

4.  Rosenzweig,  H.  "Technological  Considerations",  from  Systems  Analysis 
and  Policy  Planning:  Applications  in  Defense.  New  York:  Elsevier 
North-Roll and,Inc7,  1977. 

5.  Bate,  Roger  R. ,  et  al.  Fundamentals  of  Astrodynamlcs.  New  York: 

Dover  Publications,  Inc.,  1971. 

6.  O'Shea,  Donald  C.,  et  al.  Introduction  to  Lasers  md  Their  Applications. 
Reading,  Massachusetts:  Add'i son-Wesl ey  Publ ishing  Company,*  1978. 

7.  Havey  ,  James  H. ,  Jr.  Operational  Sciences  Department.  Lecture  materials 
distributed  in  PH  5.83,  Advanced  Weapons  Technology.  School  of  Engineer¬ 
ing,  Air  Force  Institute  of  Technology,  Wright- Patterson  Air  Force  Base, 
Ohio,  1980. 

8.  Handbook  of  Lasers  with  Selected  Data  on  Optical  Technology,  edited  by 
Pressley ,  Robert  J.  Cleveland,  Ohio":  Chemical  Rubber  Co. 

9.  Torvik, Peter  J.  Thermal  Response  Calculations  and  Their  Role  in  the 
Design  of  Experiments.  AFIT  TR  73-6.'  Air  Force  Institute  of  Technology, 
Wright-Patterson  Air  Force  Base,  Ohio,  1973. 

10,  Keeney,  Ralph  L.  and  Raiffa,  Howard.  Decisions  with  Multiple  Objectives: 
Preferences  and  Value  Tradeoffs.  New  York:  John  Wiley  &  Sons,  1976. 

11.  Haley,  Charles  and  Ghelber,  Craig.  A  Methodology  for  Validation  of 
Complex  Multi -Variable  Military  Simulation  Models.  Unpublished  Master's 
thesis.  Air  Force  Institute  of  Technology,  Wright-Patterson  Air  Force 
Base,  Ohio,  1980. 

12,  Naylor,  Thomas  H. ,  and  Finger,  J.  M.  "Verification  of  Computer  Simulation 
Models",  Management  Science,  14  (2):  B-92  to  B-101  (October  1967). 

13.  Schlesinger,  $.,  et  al.  "Developing  Standard  Procedures  for  Simulation 
Validation  and  Verification",  Proceedings  of  the  1974  Summer  Computer 
Simulation  Conference,  927-933.  La  Jolla,  California:  The~Simulation 
Counci 1 s ,  Inc. ,  1974. 


lt)9 


14.  Tytula,  Thomas  P.  A  Method  for  Validating  Missile  System  Simulation 
Models.  Technical  Report  E-Ts^Tl  for  U.  s7  Army  Missile  Research 
and  Development  Command.  Redstone  Arsenal,  Alabama:  USAMRDC,  1978. 
(AD  A055  689). 

15.  Mendenhall,  William,  and  Schaeffer,  Richard  L.  Mathematical  Statls- 
tics  with  Applications.  North  Scituate,  Massachusetts:  Dunburz 
TrMsTTmT"  — 

16.  Wilde,  Douglass  0.  Optimum  Seeking  Methods.  Englewood  Cliffs,  New 
Jersey:  Prentice-Hall,  Inc.,  1964. 

17.  Kleijnen,  Jack  P.  Statistical  Techniques  in  Simulation.  New  York: 
Marcel  Dekker,  Inc.  ,T$?5. 

18.  Pritsger,  A.  Alan  B.,  and  Pegden,  Claude  Dennis.  Introduction  to 
Simulation  and  SLAM.  West  Lafayette:  Systems  Publishing  Corporation 


19.  Shannon,  Robert  E.  Systems  Simulation--  The  Art  and  Science.  Engle 
wood  Cliffs,  New  Jersey:  Prentice-Hal  1 ,  Inc.  ,T575. 

20.  Peckham,  L.  N. ,  and  Davis,  R.  W.  "A  simplified  Propagation  Model  for 
Laser  System  Studies".  Technical  Report  Number  AFWL-TR- 72-95,  Air 
Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  1973. 

21.  Hornbeck,  Robert  W.  Numerical  Methods.  New  York,  Quantum  Publi¬ 
cations,  Inc.,  1975. 

22.  Huelsman,  Charles  B.,  Ill,  and  Barr,  David  R.  "Topics  in  Mathematics 
for  Operations  Research".  Unpublished  text  for  M.A.5.33,  Advanced 
Calculus  for  Operations  Research,  Air  Force  Institute  of  Technology, 
Wright- Patterson  Air  Force  Base,  Ohio,  1979. 


as 


A.  Date  Reduction  for.  Laser  Effects 

Representative  satellite  vulnerability  data  was  supplied  by  the 
Space  Systems  Division  of  the  Air  Force  Systems  Command's  Foreign 
Technology  Division.  Although  the  data  used  for  the  four  HELBASE  tar¬ 
get  types  is  not  "real"  data  for  any  particular  satellite  system  in  the 
real  world,  it  is  representative  of  the  "real"  data. 

This  data  was  supplied  in  the  form  shown  in  figure  22  .  The 
fluence  (in  joules  per  square  centimeter)  induced  into  the  target  is 
measured  on  the  ordinate.  The  time  (in  seconds)  that  it  takes  the 
fluence  to  be  induced  is  measured  on  the  abscissa.  A  point  on  the 
graph  represents  "failure"  of  a  satellite  subsystem  at  the  fluence  and 
time  shown. 

Since  the  information  needed  was  "time  to  satellite  negation"  at 
some  laser  intensity  (watts  per  square  centimeter),  the  data  was  trans¬ 
formed.  The  transformation  function  was 

1  ”  £  ,  where 

I  =  Laser  intensity,  (watts  per  square  centimeter) 

F  =  Fluence,  (joules  per  square  centimeter) 

t  =  Irradiation  time  (seconds). 

"I"  was  then  placed  on  the  ordinate,  with  "V  again  on  the  abscissa  (see 
figure  22  ). 

In  order  to  accomplish  this  transformation,  several  points  were  taken 
from  supplied  data  (see  table VII 1) , and  transformed.  The  problem  then  be¬ 
came  one  of  approximating  the  transformed  data  by  a  functional  represent- 
on. 
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Table  Vm 


Representative  Data  Transformation 


REPRESENTATIVE  DATA  TRANSFORMED  DATA  CONSTANTS 


TIME 

FLUENCE 

a  b  c 

TARGET  TYPE  #1 

.000.52  -.25  38 

10 

550 

55 

20 

600 

30 

50 

850 

13 

100 

1400 

14 

TARGET  TYPE  #2 

.02  -5.6  420 

10 

4200 

420 

20 

4300 

220 

50 

5000 

100 

100 

6200 

62 

TARGET  TYPE  #3 

.00012  -.7  73 

5 

500 

100 

10 

550 

55 

20 

560 

28 

50 

650 

13 

100 

780 

8 

TARGET  TYPE  #4 

.00052  -.25  38 

2 

350 

175 

5 

360 

72 

10 

370 

37 

20 

390 

20 

50 

460 

9 

100 

640 

6 
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TIME  (LOG  SCALE) 
ORIGINAL  DATA  FORM 


TIME  (LINEAR  SCALE) 
TRANSFORMED  DATA  FORM 


Fig  22. 


Representative  Data  Form 


In  all  four  cases, a  least  squares  quadratic  polynomial  was  fit 
to  the  transformed  data  (Rafs2l ; 1 22-125 •  and22 : 1 35 ) .  TablelX  also 
shows  the  quadratic  constants  generated  by  the  least  squares  quadratic 
fit. 

The  solution  for  the  negation  time  is  given  by: 

t  c  -  fe  -  [  b2  *  Hac  t  ‘iad  ^ 

where 

t  s  Irradiation  time  required  for  target  negation, 
a,  b,  c  are  quadratic  constants, 
d  is  a  constant  related  to  the  data  base, 

Ij  is  the  intensity  of  the  beam  at  the  target,  and 
Aj  is  the  wavelength  of  the  beam. 

Feasible  solutions  of  this  function  are  limited  to  times  from  1 /10th 
of  a  millisecond  to  100  seconds.  Therefore,  values  of  Ij  and  Aj  which 
would  yield  a  t  of  less  than  .0001  seconds  or  more  than  100  seconds  are 
not  accepted  by  submodel  EFFECT  (see  Appendix  C).  If  the  calculated 
negation  time  would  have  been  less  than  .0001  seconds,  it  is  assigned 
the  value  of  .0001  seconds.  If  the  calculated  negation  would  have  been 
more  than  100  seconds,  it  is  assigned  the  value  of  9,999,999  seconds. 
This  large  value  will  insure  a  target  sighting  efficiency  value  of  zero 
in  submodel  SIGHT. 


115 


B.  Recommendations 


The  recommendations  offered  here  are  divided  into  sections  describ¬ 
ing  the  general  area  of  improvement. 

Selection  of  Orbital  Weapons 

1.  Inclusion  of  elliptical  orbits.  When  only  circular  orbits  are  consider¬ 
ed  for  space  weapons,  two  of  the  six  orbital  elements  are  trivialized.  That 
is,  the  agreement  of  perigee  is  meaningless  and  the  eccentricity  is  zero. 

If  we  wish  to  consider  elliptical  orbits  in  the  generalized  optimization 
scheme,  we  therefore  pick  up  these  two  additional  degrees  of  freedom.  When 
we  have  825  orbital  possibilities  for  a  given  circular  orbit  radius,  we 
could  have  267,300  possibilities  for  a  given  elliptical  orbit  semi-major 
axis.  This  is  derived  from  a  choice  of  eccentricity  from  1.  to  .9  and  a 
5°  incremental  choice  of  agremeent  of  perigee.  (825  x  9  x  36  =  267,300). 
Due  to  the  large  number  of  orbital  possibilities,  even  searching  through 
them  for  a  "most  efficient  path"  would  be  expensive. 

Therefore,  a  random  search  technique  over  the  267,300  possibilities 
could  be  used  to  find  at  least  one  in  the  top  5%  with  a  probability  of  .9. 
This  would  require  exactly  12000  samples  of  the  267,300  possibilities. 

We  could  feel  fairly  confident  that  the  best  orbit  thus  found  would  be  near 
the  global  maximum,  and  a  multi-dimensional  search  for  the  maximum  could  be 
initiated  with  the  best  orbit  from  the  random  search  as  the  starting  point. 

The  point  is  that  extension  of  the  model  to  include  elliptical 
weapon  orbits  is  very  feasible. 

2.  Exclude  a  range  of  orbital  inclinations  from  consideration  for  a  given 
space  weapon  type.  A  feasible  scenario  for  this  option  would  involve  the 
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inability  of  available  launch  vehicles  to  place  the  weapons  into 
highly  inclined  orbits.  A  restriction  could  be  placed  on  the  SUMEFF 
matrix  and  the  search  over  the  values  by  OPTIMIZE. 

4.  Allow  orbital  "reflectors  Allowing  the  pre-hostility  placement 
of  orbital  reflectors  for  the  purpose  of  reflecting  laser  beams  to 
target  satellites  has  been  proposed.  The  submodels  SIGHT,  OPTIMIZE, 
and  BATTLE  could  be  expanded  to  include  this  interesting  concept. 

5.  Expansion  of  the  "probability  of  failure"  input.  The  simple 
probability  of  weapon  failure  (not  used  in  the  current  model)  input 
could  be  expanded  to  a  specific  reliability  function  and  failure  rate 
for  each  weapon  type.  This  would  aid  significantly  in  cost/benefit 
analyses  concerning  weapon  reliability. 

6.  Limit  the  number  of  weapons  for  a  specific  weapon  type.  Justifi¬ 
cation  for  this  option  could  be  limits  on  command  and  control  capabili¬ 
ties,  political  constraints,  or  production  limitations. 

Battle  Tactics 

Capability  to  commit  weapons  at  a  discrete  time.  In  the  case 
of'hidden"  weapons,  the  tactic  of  delaying  weapon  use  could  be  evaluated. 
This  would  only  require  the  addition  of  another  column  tc  he  TIME 
matrix  used  in  BATTLE. 

7.  Add  weapon  use  priorities.  A  possible  tactic  may  be  to  use  weapon 
X  only  against  target  types  A,  B,  and  C.  Additionally,  weapon  X  could 
be  constrained  to  be  used  first  against  target  type  B,  secondly  against 
A,  and  lastly  against  C. 

8.  Include  more  than  one  weapon  in  a  battle  station.  Common  system 
sharing  could  justify  this  option  in  many  cases. 


117 


9.  Define  a  probability  of  kill  (Pk).  A  Pk  based  on  pointing  and 
tracking  accuracies  and/or  target  maneuvering  could  be  developed  (see 
below). 

10.  Allow  target  to  maneuver  to  avoid  irradiation. 

11.  Allow  different  negation  times  for  each  target  type.  It  seems 
perfectly  feasible  that  a  user  may  desire  one  target  type  to  be  negated 
before  others.  For  example,  it  may  be  desirable  to  insure  negation  of 
target  type  A  at  the  two  hour  point,  and  type  B  at  the  ten  hour  point. 
"Mission  accomplishment"  would  require  that  both  of  these  conditions  be 
met  (or  bettered). 

12.  Weapon  slew  rate  restriction.  Calculation  of  required  weapon  slew 
rate  (cross-track  angular  velocity)  could  be  added  to  submodel  SIGHT. 
Exceeding  an  input  maximum  would  cause  "no  sighting"  for  that  encounter. 

Cost  Considerations 

13.  Revise  basic  "sighting  efficiency"  attribute.  An  attribute  to 
measure  the  "most  efficient"  weapon  path  could  be  developed  in  terms  of 
"cost  efficiency".  The  "most  efficient"  path  would  then  be  the  most 

1  ost  efficient"  path. 

As  discussed  by  several  authors,  "costing  out"  of  alternatives  appear¬ 
ed  to  be  one  of  the  few  viable  correction  factors  available  (Refs  3  : 1 25-1 27 , 
4:57-59).  It  is  plausible  that  the  "cost  per  joule"  of  an  HEL  beam  from 
a  ground-based  weapon  would  tend  to  be  much  less  than  the  "cost-per-joule" 
of  a  beam  from  a  satellite  weapon.  The  correction  factor  would  then  be 
to  "cost  out"  the  total  energy  required  to  negate  a  target.  The  relative 
energy  unit  costs  would  be  functions  of  the  HEL  energy  technology  of  each 
weapon  type. 
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So  the  attribute  to  measure  the  "most  efficient"  placement  of  weapons 
would  be  cost  effectiveness  (actually  the  reciprocal  of  cost).  Relative 
costs  could  be  assigned  to  ground-,  air-,  and  space-based  weapon  energy. 
Relative  costs  of  output  energy  must  be  functions  of  total  weapon  life 
cycle  costs  over  equal  weapon  lifetimes. 

14.  Add  fixed  cost  of  using  the  first  weapon  of  a  weapon  type.  For  some 
problems,  it  may  be  apropos  to  consider  the  relevant  fixed  costs  (research, 
development,  testing)  of  adding  the  first  weapon  of  a  previously  unused 
weapon  type  to  the  weapon  system. 
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